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Abstract 

I review works on Fe-based superconductors which depart from a metal with well defined 
Fermi surfaces and Fermi liquid-type quasiparticles. I consider normal state instabilities - SDW 
magnetism and nematic order, and superconductivity, all three as the consequences of the insta¬ 
bility of a Fermi surface due to interactions between low-energy fermionic quasiparticles. This 
approach assumes that renormalizations coming from fermions from high energies, of order band¬ 
width, modify but do not destroy Fermi liquid behavior in the normal state and can be absorbed 
into the effective low-energy model of interacting fermions located near hole and electron-type 
Fermi surfaces. I argue that the interactions between these fermions are responsible for (i) a 
stripe-type SDW magnetic order (and, in some special cases, a checkerboard order ) , (ii) a 
pre-emptive nematic-type instability, in which magnetic fluctuations break C 4 lattice rotational 
symmetry down to C 2 , but magnetic order does not yet develop, and (hi) a superconductivity, 
which competes with these two orders. The experimental data on superconductivity show very 
rich behavior with potentially different symmetry of a superconducting state even for different 
compositions of the same material. I argue that, despite all this, the physics of superconduc¬ 
tivity in the itinerant scenario for Fe-based materials is governed by a single underlying pairing 
mechanism. 
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I. INTRODUCTION 


The discovery of superconductivity in Fe-based pnictides [1] (Fe-based compounds with 
elements from the 5th group: N, P, As, Sb, Bi) was, arguably, among the most signihcant 
breakthroughs in condensed matter physics during the past decade. A lot of efforts by 
the condensed-matter community have been devoted in the few years after the discovery 
to understand normal state properties of these materials, the pairing mechanism, and the 
symmetry and the structure of the pairing gap. 

The family of Fe-based superconductors (FeSCs) is already quite large and keeps 
growing. It includes various Fe-pnictides such as fill systems RFeAsO {R =rare earth 
element) [1-4], 122 systems XFe 2 As 2 (X=alkahne earth metals) [5-7], 111 systems like 
LiFeAs [8], and also Fe-chalcogenides (Fe-based compounds with elements from the 16th 
group: S, Se, Te) such as FeTei_a;Sej, [9] and Aj,Fe 2 -ySe 2 (A = K,Rb,Cs) [10,11]. 

Superconductivity (SC) in FeSCs emerges upon either hole or electron doping (see 
Fig. 1), but can also be induced by pressure or by isovalent replacement of one pnictide 
element by another, e.g.. As by P (Ref. [12]). In some systems, like LiFeAs [8], LiFeP 
[13] and LaFePO [14], SC emerges already at zero doping, instead of a magnetic order. 

Parent compounds of nearly all FeSCs are metals, in distinction to cuprate supercon¬ 
ductors for which all parent compounds are Mott insulators. Still, in similarity with the 
cuprates, in most cases these parent compounds are antiferromagnetically ordered [15]. 
Some researchers [16-18] used this analogy to argue that FeSCs are at short distance from 
Mott transition, and at least some elements of Mott physics must be included into the de¬ 
scription of these systems. A rather similar point of view is [18] that fermionic excitations 
in FeSCs display both localized and itinerant properties and the interplay between the 
two depends on the type of the orbital (one set of ideas of this kind lead to the notion of 
"orbital selective Mott transition on FeSCs [17,18]). An alternative point of view, which 
I will present in this review, is that low-energy properties of most of FeSCs can be fully 
captured in a itinerant approach, without invoking Mott physics. 

In itinerant approach, electrons, which carry magnetic moments, travel relatively freely 
from site to site. The magnetic order of such electrons is often termed as a spin-density- 
wave (SDW), by analogy with e.g., antiferromagnetic Cr, rather than "Heisenberg anti¬ 
ferromagnetism" - the latter term is reserved for systems in which electrons are "nailed 
down" to particular lattice sites by very strong Coulomb repulsion. From experimental 
perspective, the majority of FeSCs display a rather small ordered moment in the normal 
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state, consistent with SDW scenario [19]. There are notable exceptions - Fe-chalcogenide 
FeTe (the parent compound of FeTei_j,Sea;, which superconduct at x around 0.5) displays 
magnetic properties consistent with the Heisenberg antiferromagnetism of localized spins 
[20]. However, the properties of this material vary quite substantially between a; = 0 
and X = 0.5, and magnetic fluctuations at a; ~ 0.5 are similar to those of other FeSCs. 
Another example where magnetism is strong and probably involves localized carriers is 
Aj,Fe 2 -ySe 2 (Ref. [10]). However, in this material, localized carriers and itinerant carriers 
are most likely phase separated, with superconductivity coming primarily from itinerant 
carriers. 

The itinerant approach to magnetism and superconductivity in FeSCs and the compar¬ 
ative analysis of Fe- and Cu-based superconductors have been reviewed in several recent 
publications [19,21-32]. This review is an attempt to summarize our current understand¬ 
ing of the phase diagram, the origin of SDW and nematic orders, the pairing mechanism 
for superconductivity, and the symmetry and the structure of the pairing gap at various 
hole and electron dopings. 

Like I said, the very idea of itinerant approach is that magnetism and superconductivity 
come from the interactions between fermionic states located very near the Fermi surfaces. 
These interactions originate from a Coulomb interaction, which is obviously a repulsive 
one. 

A repulsive interaction between itinerant carriers is well known to lead to Stoner-type 
magnetic instability, and the presence of the SDW-ordered phase on the phase diagram 
of FeSCs should not come as a surprise. Less obvious issue is what kind of magnetism is 
present in FeSCs. Experiments show that most of undoped and weakly doped Fe-pnictides 
display the stripe spin-density wave order at T < Tgdw) with ordering vectors ( 0 , 7 r) or 
(tt, 0) in the 1-Fe Brillouin zone (IFeBZ), Ref. [33,34] (see Fig. 2). Such an order not only 
breaks 0(3) spin symmetry, but also breaks lattice rotational symmetry from O 4 down 
to O 2 (the stripes run either along X or along Y direction). Stripe, order, however, does 
not emerge in all cases. Neutron scattering data on more heavily doped Bai_a;Naa;Fe 2 As 2 
(Ref. [ 35]) and on Ba(Fei_a,Mn 2 .) 2 As 2 (Ref. [36]) show that the magnetic order there 
does not break O 4 symmetry (examples are shown in Fig. 2). I will argue that both types 
of magnetic order (the one which breaks 6*4 symmetry and the one which doesn’t) emerge 
in the itinerant scenario for FeSCs. 

Another interesting aspect of the normal state phase diagram is that in weakly doped 
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Figure 1: Schematic phase diagram of Fe-based pnictides upon hole or electron doping. In the 
shaded region, superconductivity and antiferromagnetism co-exist. Not all details/phases are 
shown. Superconductivity can be initiated not only by doping but also by pressure and/or 
isovalent replacement of one pnictide element by another [12]. Nematic phase at T > T^r 
is subject of debates. Superconductors at large doping are KFe 2 As 2 for hole doping [86,89] 
and A 2 ;Fe 2 -ySe 2 (A = K, Rb, Cs) for electron doping [10,11]. Whether superconductivity in 
pnictides exists at all intermediate dopings is not clear yet. From Ref. [29]. 

Fe-pnictides, the stripe SDW order is often preceded by a “nematic” phase with broken 
C 4 tetragonal symmetry but unbroken 0(3) spin rotational symmetry. The emergence 
of such a phase is not only manifested by a tetragonal to orthorhombic transition at 
Tn > Tsdw, but also by the onset of signihcant anisotropies in several quantities [37], such 
as dc resistivity [38,39], optical conductivity [40,41], local density of states [42], orbital 
occupancy [43], susceptibility [44], and the vortex core in the mixed superconducting state 
[45]. The fact that the SDW and structural transition lines follow each other across all the 
phase diagrams of 1111 and 122 materials, even inside the superconducting dome [46,47], 
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Figure 2: (color online) Various SDW spin configurations described by For 

a perfect nesting only A^ + A^ is fixed. Panel (a) - Ai = 0, panel (b) - A 2 = 0, panel (c) - 
Ai A A 2 , and panel (d) - Ai = A 2 . From [120] 























prompted researchers to propose that SDW and nematic orders are intimately connected. 
The interplay between magnetic and structural transitions in FeSCs is also quite rich: 
while in 1111 materials the two transitions are second-order and split (T„ > Tcdw), in most 
of the 122 materials they seem to occur simultaneously or near-simultaneously at small 
dopings, but clearly split above some critical doping - a; ~ 0.022 in Ba (Fei_ 3 ;Coa ;)2 As 2 , 
see [48,49], and x ~ 0.039 in Ca (Fei_ 2 ;Coa ,)2 As 2 , see [50]. 

For superconductivity, the central issue is what causes the attraction between fermions. 
The BCS theory of superconductivity attribute the attraction between fermions to the 
underlying interaction between electrons and phonons [51] (the two electrons effectively 
interact with each other by emitting and absorbing the same phonon which then serves as 
a glue which binds electrons into pairs). Electron-phonon mechanism has been successfully 
applied to explain SC in a large variety of materials, from Hg and Al to recently discovered 
and extensively studied MgB 2 with the transition temperature = 39K [52]. However, 
for FeSCs, early hrst-principle study of superconductivity due to electron-phonon inter¬ 
action placed Tc at around IK, much smaller that the actual Tc in most of FeSCs. This 
leaves an electron-electron interaction as the more likely source of the pairing. 

Pairing due to electron-electron interaction has been discussed even before high Tc era, 
most notably in connection with superfluidity in ^He [53,54], but became the mainstream 
after the discovery of SC in the cuprates [55]. This discovery signaled the beginning of 
the new era of “high-temperature superconductivity” to which FeSCs added a new avenue 
with quite high traffic over the last hve years. 

A possibility to get superconductivity from nominally repulsive electron-electron in¬ 
teraction is based on two fundamental principles. First, in isotropic systems the analysis 
of superconductivity factorizes [56] between pairing channels with different angular mo¬ 
menta / = 0,1, 2, 3, etc [in spatially isotropic systems I = 0 component is called s—wave, 
I = 1 component is called p—wave, I = 2 component is called d—wave, and so on]. If 
just one component with some I is attractive, the system undergoes a SC transition at 
some temperature T = Tc. Second, the screened Coulomb interaction U{r) is constant 
and repulsive at short distances, but oscillates at large distances and may develop an 
attractive component at some 1. Kohn and Luttinger (KL) have explicitly proven back 
in 1965 (Ref. [57]) that the combination of these two effects necessary leads to a pairing 
instability, at least at large odd /, no matter what the form of U (r) is. 

In lattice systems, angular momentum is no longer a good quantum number, and the 
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equation for only factorizes between different irreducible representations of the lattice 
space group. In tetragonal systems, which include both cuprates and FeSCs , there are 
four one-dimensional irreducible representations Aig, Big, i? 2 g, and A 2 g and one two- 
dimensional representation E2g. Each representation has inhnite set of eigenfunctions. 
The eigenfnnctions from Aig are invariant nnder symmetry transformations in a tetragonal 
lattice: x —)■ —x, y —)■ —y, x ^ y, the eigenfnnctions from Big change sign nnder x ^ y, 
and so on. If a snpercondncting gap has Aig symmetry, it is often called s—wave because 
the hrst eigenfunction from Aig gronp is jnst a constant in momentnm space (a <5—fnnction 
in real space). If the gap has Big or B2g symmetry, it is called d—wave {dx2-y2 or dxy, 
respectably), because in momentum space the leading eigenfunctions in Big and B2g are 
cos kx — cos ky and sin kx sin ky, respectively, and these two rednce to / = 2 eigenfnnctions 
cos 2d and sin 2d in the isotropic limit. 

In the cnprates, the snpercondncting gap has been proved experimentally to have Big 
symmetry [58]. Such a gap appears quite naturally in the doping range where the cnprates 
are metals, because KL-type consideration shows that Big interaction becomes attractive 
if the fnlly dressed repnlsive interaction between fermions near different corners of the 
Brillonin zone (the one at momentnm transfer near (tt, tt)) exceeds the repnlsion at small 
momentnm transfer. The enhancement of (tt, tt) interaction is a snre thing if the system 
displays strong antiferromagnetic spin flnctnations (see Fig.3). That Big gap is selected 
is not a snrprise because such gap A{k) oc coskx — cosky changes sign not only nnder 
kx —t ky bnt also between k and k' = k -|- Q where Q = (tt, tt). This sign change is the 
crncial element for any electronic mechanism of snpercondnctivity because one needs to 
extract an attractive component from repnlsive screened Conlomb interaction. 

In FeSCs, magnetism and snpercondnctivity are also close neighbors on the phase 
diagram, and it has been proposed [59,60] at the very beginning of the Fe era that the 
pairing mechanism in FeSCs is also a spin-flnctnation exchange. However, the geometry 
of low-energy states in FeSCs and in the cnprates is different, and in most FeSCs the 
momentnm Q connects low-energy fermionic states near the center and the corner of the 
Briilonin zone (see Fig.3). A simple experimentation with trigonometry then tell us that 
the SC order parameter (the gap) A(/c) must be symmetric with respect to kx —)■ ky 
and kx —)■ —kx, bnt still must change sign nnder k —>■ k -T Q. Snch order parameter 
belongs to Aig representation, but it only has contribntions from a particular subset of 
Aig states with the form coskx + cosky, cosSkx + cos(3/cy), etc, which all change sign 
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Figure 3: A comparison of the pairing state from spin-fluctuation exchange in cuprate SCs and 
in FeSCs. In the cuprates (left panel) the FS is large, and antiferromagnetic Q = (tt, tt) connects 
points on the same FS. Because spin-mediated interaction is positive (repulsive), the gap must 
change sign between FS points separated by Q. As the consequences, the gap changes sign 
twice along the FS. This implies a d—wave gap symmetry. In FeSCs (left panel) scattering by Q 
moves fermions from one FS to the other. In this situation, the gap must change sign between 
different FS, but to first approximation remains a constant on a given FS. By symmetry, such 
a gap is an s—wave gap. It is called s~^ because it changes sign between different FSs 

under k —k-|- Q. An order parameter with such symmetry is called an extended s—wave 
or, in shorter notations, 

Majority of researches do believe that in weakly/moderately doped FeSCs the gap 
does have symmetry. However, numerous studies of superconductivity in FeSCs over 
the last hve years demonstrated that the physics of the pairing is more involved than it 
was originally thought because of multi-orbital/multi-band nature of low-energy fermionic 
excitations in FeSCs. It turns out that both the symmetry and the structure of the pairing 





































gap result from rather non-trivial interplay between spin-fluctuation exchange, intraband 
Coulomb repulsion, and momentum structure of the interactions. In particular, an gap 
can be with or without nodes, depending on the orbital content of low-energy excitations. 
Besides, the structure of low-energy spin fluctuations evolves with doping, and the same 
spin-fluctuation mechanism that gives rise to s’*"” gap at small/moderate doping in a 
particular material can give rise to a d—wave gap at strong hole or electron doping. 

There is more uncertainly on the theory side. In addition to spin fluctuations, FeSCs 
also possess charge fluctuations whose strength is the subject of debates. There are 
proposals [61,62] that in multi-orbital FeSCs charge fluctuations are strongly enhanced 
because the system is reasonably close to a transition into a state with an orbital order 
- a spontaneous symmetry breaking between the occupation of different orbitals). (A 
counter-argument is that orbital order does not develop on its own but is induced by 
a magnetic order [63,64]). If charge fluctuations are relevant, one should consider, in 
addition to spin-mediated pairing interaction, also the pairing interaction mediated by 
charge fluctuations. The last interaction gives rise to a conventional, sign-preserving 
s—wave pairing [61]. A "p-wave" gap scenario (a gap belonging to E 2 g representation) 
has also been put forward [65]. 

From experimental side, s-wave gap symmetry is consistent with ARPES data on 
moderately doped Fe 2 As 2 and BaFe 2 (Asi_a;Pa;) 2 , which detected only a small 

variation of the gap along the FSs centered at (0, 0) (Ref. [66]), and with the evolution of 
the tunneling data in a magnetic field [67]. However, other data on these and other FeSCs, 
which measure contributions from all FSs, including the FSs for which ARPES data are 
not available at the moment, were interpreted as evidence either for the full gap [68-71], 
or that the gap has either accidental nodes [72,73] or deep minima [74-76]. As additional 
level of complexity, supercondnctivity was also discovered in materials which only contain 
hole pockets, like hole-doped KFe 2 As 2 , or only electron pockets, like A 2 ;Fe 2 - 3 /Se 2 . For 
these materials, the argument for s’*"" snperconductivity, driven by magnetically-enhanced 
interaction between fermions near hole and electron pockets, is no longer applicable, yet 
both classes of materials have hnite T^, which is around 3K for KFe 2 As 2 and as high as 30K 
for Aa,Fe 2 -j,Se 2 (Refs. [77]). For KFe 2 As 2 , Various experimental probes [89] indicate the 
presence of gap nodes. Laser ARPES data [78] were interpreted as evidence for s—wave 
with nodes, while thermal condnctivity data have been interpreted as evidence for both, 
d—wave and s—wave orders (Refs.[79] and [80], respectively). For A 2 ;Fe 2 -j;Se 2 , ARPES 
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results were interpreted as evidence for s—wave (Ref. [11]), however neutron scattering 
experiments [81] detected a resonance peak which most naturally can be interpreted as 
evidence for d—wave [82] (see, however, [83]). 

In this paper, I argue that all these seemingly very different gap structures actually 
follow quite naturally from the same underlying physics idea that FeSCs can be treated as 
moderately interacting itinerant fermionic systems with multiple FS sheets and effective 
four-fermion intra-band and inter-band interactions in the band basis. I introduce the 
effective low-energy model with small numbers of input parameters [84] and use it to 
study the doping evolution of the pairing in hole and electron-doped FeSCs. I argue that 
various approaches based on underlying microscopic models in the orbital basis reduce to 
this model at low energies. 

The paper is organized as follows. In Sec. II I discuss general aspects of the band 
structure of FeSCs which contain hole and electron pockets. In Sec. Ill I present a 
generic discussion of what is needed for SDW order and superconductivity and how mag¬ 
netic fluctuations help superconductivity to develop. In Sec. IV I briefly review parquet 
renormalization group approach to FeSCs. This approach treats magnetism and supercon¬ 
ductivity on equal footing. I argue that, depending on input parameters and/or doping, 
the system first becomes either SDW magnet or a superconductor. In Sec.VI I review 
itinerant approach to magnetism. I show that for most (but not all) dopings a SDW order 
below Tsdw spontaneously breaks C 4 lattice symmetry in addition to 0(3) symmetry of 
rotations in spin space. I then review works on a pre-emptive spin-nematic instability at 
Tn > Tsdw, when the system spontaneously breaks O 4 symmetry down to O 2 , but spin- 
rotational symmetry remains unbroken down to a smaller Tgdw- In Sec.VII I review an 
itinerant approach to superconductivity. I hrst present a generic symmetry considera¬ 
tion of a gap structure in a multi-band superconductor and show that a “conventional 
wisdom” that an s-wave gap is nodeless along the FSs, d-wave gap has four nodes, etc, 
has only limited applicability in multi-band superconductors, and there are cases when 
the gap with nodes has an s—wave symmetry, and the gap without nodes has a d—wave 
symmetry. I then discuss the interplay between intra-band and inter-band interactions, 
for realistic multi-pocket models for FeSCs and set the conditions for an attraction in an 
s—wave or a d—wave channel. I consider 5-orbital model with local interactions, con¬ 
vert it into a band basis, and show the structure of the superconducting gap. I use the 
combination of RPA and leading angular harmonic approximation to analyze the pairing 
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Figure 4: Crystallographic structures of various families of iron-based superconductors. From 
Ref. [21]. 

in s— and d—wave channels at different dopings. I show that, depending on parameters 
and doping, magnetically-mediated pairing leads to an superconductivity with either 
a near constant gap along the FSs, or gaps with deep minima, or even with the nodes. I 
briefly review the experimental situation in Sec. VIII and present concluding remarks in 
Sec. IX. 


II. THE ELECTRONIC STRUCTURE OF FESCS. 

The crystallographic structures of various families of iron-based superconductors is 
shown in Fig. 4. All FeSCs contain planes made of Fe atoms, and pnictogen/chalcogene 
atoms are staggered in a checkerboard order above and below the iron planes. In fill 
system this order repeats itself from one Fe plane to the other, while for 122-type systems, 
it flips sign between neighboring planes. 

The electronic structures of FeSCs at low energies are rather well established by 
ARPES [85] and quantum oscillation measurements [87]. In weakly and moderately 
electron-doped materials, like BaFei_a;Coa,Fe 2 As 2 the FS contains several quasi-2D warped 
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Figure 5: The electronic structure of FeSCs. In weakly and moderately electron-doped materials 
(left panel) the FS consists of quasi-2D warped cylinders centered at (0, 0) and (vr, vr) in a 2D 
cross-section. The ones near (0, 0) are hole pockets (filled states are outside cylinders), the ones 
near (vr, tt) are electron pockets (filled states are inside cylinders) There also exists a quasi-3D 
hole pocket near kz = vr. In hole-doped FeSCs the electronic structure is very similar, but 3D 
hole pocket becomes quasi-2D warped hole cylinder. From Ref. [28]. 

cylinders centered at fc = (0, 0) and k = (vr, vr) in a 2D cross-section, and may also contain 
a quasi-3D pocket near kz = ir (Fig.5). The fermionic dispersion is electron-like near the 
FSs at (vr, vr) (hlled states are inside a FS) and hole-like near the FSs centered at (0, 0) 
(hlled states are outside a FS). In heavily electron-doped FeSCs, like A 3 ;Fei_j^Se 2 (A = K, 
Rb, Cs), only electron pockets remain, according to recent ARPES studies. [10] In weakly 
and moderately hole-doped FeSCs, like Bai_a:Ka,Fe 2 As 2 , the electronic structure is similar 
to that at moderate electron doping, however the spherical FS becomes the third quasi 
2D hole FS centered at (2vr, 0) = (0,0). In addition, new low-energy hole states likely 
appear around (vr, vr) and squeeze electron pockets [88]. At strong hole doping, electron 
FSs disappear and only only hole FSs are present [86] These electronic structures agree 
well with hrst-principle calculations [22,91,92], which is another argument to treat FeSCs 
as itinerant fermionic systems. 

The measured FS reflects the actual crystal structure of FeSCs in which there are two 
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Figure 6: Upper panel: 3D electronic structure of LaOFeAs (left) and its 2D cross-section (left). 
In only Fe states are considered, an elementary cell contains one Fe atom (green). The actual 
unit cell (blue) contains two Fe atoms because of two non-equivalent positions of a pnictide 
above and below the Fe plane. Lower panel - the location of hole and electron FSs in a 2D cross 
section in the folded BZ (two Fe/cell, right) and in the unfolded BZ (one Fe/cell, left). From 
Refs. [191], [192](b) and [193]b. 

non-equivalent positions of a pnictide above and below an Fe plane, and, as a result, 
there are two Fe atoms in the unit cell (this actual 2Fe BZ is called "folded BZ"). From 
theory perspective, it would be easier to work in the BZ which contains only one Fe 
atom in the unit cell (this theoretical IFe BZ is called "unfolded BZ"). I illustrate the 
difference between folded and unfolded BZ in Fig.6. In general, only folded BZ is physically 
meaningful. However, if by some reason a potential from a pnictogen (or chalcogen) can 
be neglected, the difference between the folded and the unfolded BZ becomes purely 
geometrical: the momenta kx and ky in the folded BZ are linear combinations of kx and 
ky in the unfolded BZ: kx = kx + ky, ky = kx — ky. In this situation, the descriptions in 
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the folded and unfolded BZ become equivalent. 

Most of the existing theory works on magnetism and on the pairing mechanism and 
the structure of the SC gap analyze the pairing problem in the unfolded BZ, where which 
two hole pockets are centered at (0,0) and one at (7r,7r), and the two electron pockets 
are at (0,7r) and (7r,0). It became increasingly clear recently that the interaction via a 
pnictogen/chalcogen and also 3D effects do play some role for the pairing, particularly 
in strongly electron-doped systems. [93,94] However, it is still very likely that the key 
aspects of the pairing in FeSCs can be understood by analyzing a pure 2D electronic 
structure with only Fe states involved. In the next three sections I assume that this is 
the case and consider a 2D model in the unfolded BZ with hole FSs near (0, 0) and (tt, tt) 
and electron FSs at (0,7r) and (tt, 0). 

III. THE LOW-ENERGY MODEL AND THE INTERPLAY BETWEEN MAG¬ 
NETISM AND SUPERCONDUCTIVITY 

For proof-of-concept I first consider a simple problem: a 2D two-pocket model with 
one hole and one electron FS, both circular and of equal sizes (see Fig.7), and momentum- 
independent four-fermion interactions. 

The free-fermion Hamiltonian is the sum of kinetic energies of holes and electrons: 

H2 = Y. ^cc{^Ck,a + effl^fk,a ( 1 ) 

k,a 

where c stands for holes, / stands for electrons, and Ecj stand for their respective dis¬ 
persions with the property that edk) = —ef{k + Q), where Q = (tt, tt) is the momentum 
vector which connects the centers of the two fermi surfaces. The density of states No is 
the same on both pockets, and the electron pocket ‘nests’ perfectly within the hole pocket 
when shifted by Q. 

There are hve different types of interactions between low-energy fermions: two intra¬ 
pocket density-density interactions, which I treat as equal, interaction between densities 
in different pockets, exchange interaction between pockets, and pair hopping term, in 
which two fermions from one pocket transform into two fermions from the other pocket. 
I show these interactions graphically in Fig 8. 
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Figure 7: FS topology for a two-pocket model for FeSCs. The two FSs are for hole-like dispersion 
(blue circle, filled states outside the FS) and electron-like dispersion (orange circle, filled states 
inside the FS). 



Figure 8 : The interactions between pockets in the two-pocket model for Fe-pnictides. Gi is 
a density-density interaction between fermions from different pockets. G 2 is an exchange in¬ 
teraction between the pockets, G 3 is a pair hopping process between the pockets, and G^ is 
a density-density interaction within the same pocket. All interactions are repulsive (positive). 
From [32,101]. 
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Figure 9: Vertices F/j/j = Fee and F/^e introduced in the 2 pocket model. Solid and dashed lines 
denote fermions from the two pockets. From [112]. 

In explicit form 

H'uit = Gl ^kiafk2cr'fk3<j'Ck4ij 

[k,a] 

^ ^ fk\a^k2cr'fk^u'^k^a 
[k,a] 

^ki,iJi^k2,iJ2fk3,(T2fk4,ai + h.C^ 

[k,(T] 

+ E 

[fc.O-] 

( 2 ) 

where [k,o-] is short for the sum over the spins and the sum over all the momenta 
constrained to ki + k 2 = + k^ modulo a reciprocal lattice vector. 

The textbook approach to analyse potential instabilities towards superconductivity 
and magnetism is to consider the appearance of the poles in the corresponding vertex 
functions. For superconductivity, we need to consider vertex functions with zero total 
incoming momentum; —kj?; pi?, —pi?); ree(kir, —kp; p^?, —Pf), where k^ and Pf 

belong to the same pocket, and r/ie(kF, —kF; Pf, —Pf), where kF and pf belong to 


Ct4 I 

/O I 


^ki,ai ^k2,cr2 *^^3F2*^^4FI 


+ C^ /) 
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Figure 10: Lowest-order terms in the ladder series for the renormalizations of the SDW 
superconducting vertices. The effective coupling in the SDW channel is Gi + G 3 . The 
couplings in the two SDW channels are G 3 + G 4 and — 6*3 + G 4 . From [32,101]. 
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different pockets (see Fig. 9). To first order in Gj, we have 


—k^; Pf, —Pf) = —Gi 
^eei^F, -kF; PF, “Pf) = -G4 

kL(kF, —kF; Pf, — Pf) = —G3 (3) 

I follow^® and introduce the vertex function with the opposite sign compared to the 
interaction potential. 

For SDW order we need to consider interactions with momentum transfer Q = (7r,7r): 
rL(kF, k),; Pf, Pf), F5)e(kF, pV; k'F, Pf), and F5)e(kF, pV; Pf, k),), where k and k' belong 
to one pocket and p and p' belong to the other pocket, and p = k + Q. To first order in 
G we have 


rL(kF, k);^; Pf, Pf) = —G3 
rL(kF,p'F;kF,PF) = -Gi 

rL(kF, P f; Pf, k),) = -G 2 (4) 

To see which combinations of different F appear in the SDW and superconducting 
channels, I add to the Hamiltonian the trial terms Aic[ a'^ai 3 fk+Q,/ 3 , and 

Aefl ^iayai 3 f-k,f 3 , dress them by the interactions, and express the fully renormalized Ai, 
Ae, and via fully renormalized vertices. The lowest-order terms in the corresponding 
series are shown in Fig. 10. One can easily make sure that the vertex which renormalizes 
Ai contains Gi + G 3 , while the vertices which renormalize A^ and A^ are made out of 
G 3 and G 4 . 


A. Ladder approximation 

To proceed further, I hrst assume that the two channels do not communicate with each 
other, i.e., the renormalization of the SDW vertex does not involve the interactions with 
zero total momentum, while the renormalization of the two superconducting vertices does 
not involve the interaction with momentum transfer Q. Mathematically, this approxima¬ 
tion implies that higher-order additions to Fig. 10 form ladder series. These series can 
be easily summed up analytically. 
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1. The SDW vertex 


For SDW vertex, summing up ladder diagrams we obtain 

Ai = Af>v.«). XpftW) = (i - n,4Q)rLt) (s) 

where 

-pfuii _ _ Gi + _ 

l-Up,{Q){G, + Gs) 

where Ilph{Q) is the particle-pole polarization bubble at momentum transfer Q. Note 
that only the combination Gi + G 3 appears in ( 6 ). The interactions G 2 and G 4 do not 
participate in the renormalization of the SDW vertex. 

I show the behavior of Uphiq) at a generic q in Fig. 12 below. At this stage, it is just 
enough to observe that Ilph{Q) is positive. Eq. ( 6 ) then shows that the full vertex in the 
SDW channel and the susceptibility Xph{Q) diverge when np/i(( 5 )(Gi-l-G 3 ) = 1. That 
the divergence occurs for a repulsive interaction {Gi -|- G 3 > 0 ) reflects the well-known 
fact that fermion-fermion repulsion does give rise to a Stoner-like magnetic instability. 


2. The supereondueting vertex 


Let’s now solve for the full Ah and Ag in the ladder approximation. A simple analysis 
shows that the two equations become 


A* = Af - (Ai.“'r£" + AfX'f) n„ 

A. = Af' - {A'prif + A<"'r(f) n„ (7) 


where Upp > 0 is the particle-particle polarization bubble at zero momentum transfer: 
(flpp = Ai"o(log IcUc/D] + where Nq is the density of states at the Fermi level and D 

is the total incoming frequency), and 


full 
^ hh 


-pfull 


full 
^ he 


_1 ( G, + G3 G,-G3 \ 

2 \1 + (G^4 + G3)npp 1 -l- (G 4 — j 

yn full 

_1 ( G4 + G3 _ G4 - G3 \ 

2 v 1 + (G4 -i- G3)npp 1 + (G4 — G3)npp j 


( 8 ) 
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The set of equations in (7) decouples into 


Xpp I ^ 

A;, + Ae = (aI ) + A®) x+, xtp = 1 + (G 4 + G 3 )Upp 

(9) 

Because IIpp > 0, the presence or absence of a pole in (i.e., potential divergence of 
Xpp) depends on the signs of G 3 + G 4 or G 4 — G 3 . If both are positive, there are no poles, 
i.e., non-superconducting state is stable. In this situation, at small Q, k, -i/n,p, 
~ — Gl))Ilpp, i.e., both vertex functions decrease (inter-pocket vertex 

decreases faster). If one (or both) combinations are negative, there are poles in the upper 
frequency half-plane and fermionic system is unstable against pairing. The condition for 
the instability is IG 3 I > G 4 . G 4 is inter-pocket interaction, and there are little doubts 
that it is repulsive, even if to get it one has to transform from orbital to band basis. G 3 
is interaction at large momentum transfer, and, in principle, it can be either positive or 
negative depending on the interplay between intra- and inter-orbital interactions. In most 
microscopic mnlti-orbital calculations, G3 tnrns ont to be positive, and I set G3 > 0 in 
the analysis (for the case G 3 < 0 see Ref.®®). 

For positive G 3 , the condition for the pairing instability is G 3 > G 4 . What kind of 
a pairing state do we get? First, both and F^g^^ do not depend on the direction 
along each of the two pockets, hence the pairing state is necessary s—wave. On the other 
hand, the pole is in F 2 , which appears with opposite sign in and F^g^^ The pole 
components of the two vertex functions then also differ in sign, which implies that the 
two-fermion pair wave fnnction changes sign between pockets. Snch an s—wave state is 
often call s’*"” to emphasize the sign change between the pockets. This wave fnnction 
mnch resembles the second wave fnnction from Aig representation: coskx + cosky. It is 
still s—wave, bnt it changes sign under k —^ k -|- (tt, tt), which is precisely what is needed 
as hole and electron FSs are separated by {71,71). I cantion, however, that the analogy 
should not be taken too far becanse the pairing wave fnnction is dehned only on the two 
FSs, and any fnnction from Aig representation which changes sign under k —^ k -|- (tt, vr) 
would work equally well. 
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B. Beyond ladder approximation 


1. How to get an attraction in the pairing channel? 

Having established the pairing symmetry, I now turn to the central issue: how to get 
an attraction in the pairing channel? Let’s start with the model with a momentum- 
independent (Hubbard) interaction in band basis. For such interaction, all Gi are equal, 
i.e, G 3 = G 4 = Gi = G. The SDW vertex still diverges when 2GYiph{Q) = 1, but Xpp = 1 
and Xpp vanishes at small This implies that, within ladder approximation, the only 
instability is a SDW. This does not holds, however, beyond the ladder approximation, as 
I now demonstrate. The consideration below follows Refs.^^’^°^. 

Kohn-Luttinger consideration 

As the hrst step away from the ladder approximation, consider how KL physics works 
in our case. By this I mean that the intra-pocket interaction G 4 and pair-hopping G 3 are 
both equal to G only if they are treated as bare interactions. In reality, each of the two 
should be considered as irreducible interaction in the pairing channel. The irreducible 
interaction is the bare interaction plus all renormalizations except for the ones in the 
particle-particle channel. KL considerations includes such renormalizations to order G^. 
Below I label irreducible pairing vertices as T^/i and T^g. 

The contributions to T^^j and T^g to order G^ are shown in Fig 11. In analytical form 
I have 

= -G 4 - {gI + Gl- 2Gi(G'i - G 2 )) np,(0), 

fL = -G 3 - 2G'3(2Gi - G2)nph{Q), (10) 

where, I remind, Q = (tt, tt). For a constant G this reduces to 

= -G(l + 2Gnpft(0)), 

fL = -G(l + 2Gn,^((3)), (11) 

One can show that the relation (8) still holds if we replace G 3 by —F°g and G 4 by —F^^^. 
Because F{g“ = I will only deal with and F{g^\ which are given by 

1 / -I- ro \ 

“ ^ 2 L-(fL + fgi)n„ ' 

1 / ro _L ro _ ro \ 

_ _ I _ ^ he ' ^ hh _ _ _ ^ hh ^ he _ j 

- 2 L-(fL + fgi)n„ 1 - (fj, - fL)n„l ^ 

( 12 ) 
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contributions with 11(0) 


contributions with n(Q) 




Figure 11: Contributions to the irreducible vertices r°^(top) and F/ie(bottom). F®^ only gets 
contributions form 11(0) while F^g gets contribution from n(Q). From [32]. 
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The condition for the pairing instability becomes |r°g| > —T^^. Comparing the two 
irreducible vertex functions, I hnd 

+ IfLl = 2G^ (n„.(Q) - n„.(0)) (13) 

i.e., the condition for the pairing is satished when Yiph{Q) > np/i(0). For a gas of fermions 
with one circular FS, Iiph{q) either stays constant or decreases with g, and the condition 
^phiQ) > np;i(0) cannot be satished. However, in our case, the two FS’s are separated 
by Q, and, moreover, one FS is of hole type, while the other is of electron type. One can 
easily verify that, in this situation, Yiph{Q) is enhanced comparable to np/i(0). I present 
the plot of Ylph{q) along = qy in Fig 12. Indeed, Ylph{Q) is much larger than np/j(0). 

We see therefore that for the renormalization of the bare interaction into an irreducible 
pairing vertex does give rise to an attraction in the Aig pairing channel. The attractive 
pairing interaction is weak and at this stage is certainly smaller than the interaction 
in the SDW channel. On the other hand, the polarization bubble 'n.ph{Q) is in general 
some constant, while the polarization bubble Tlpp diverges logarithmically when the total 
frequency O vanishes. 

Before I proceed, a comment. Because we deal with fermions with circular FSs lo¬ 
cated near particular /c—points, polarization operators at small momentum transfer and 
momentum transfer Q = (tt, tt) can be approximated by constants. Then the irreducible 
vertex function has only an s—wave {Aig) harmonic, like the bare interaction, i.e. KL 
renormalization does not generate interactions in other channels. Treating pockets as 
circular is indeed an approximation, because for square lattice the only true requirement 
is that each FS is symmetric with respect to rotations by multiples of 7 r /2 symmetry). 
For small pocket sizes, deviations from circular forms are small, but nevertheless are gen¬ 
erally hnite. If we include this effect, we hnd that the KL effect does generate interactions 
in other channels (i?ig,i? 2 g, and A 2 g), which may be attractive, and also leads to more 
complex structure of the pair wave function in s^~ channel, which now acquires angular 
dependence along hole and electron pockets, consistent with C 4 symmetry®®’®^ 

The Hubbard limit of a constant G is a somewhat artihcial case, however. The actual 
bare interactions Gi have to be extracted from the multi-orbital model and do depend on 
momentum transfer. In this situation G 4 — G 3 is generally non-zero already before KL 
renormalization. It is natural to expect that the bare interaction is a decreasing function 
of momenta, in which case G 4 , which is the interaction at small momentum transfer, is 
larger than the interaction G 3 at momentum transfer near Q. Then the KL term has to 
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compete with the first-order repulsion. As long as G'n.ph{Q) is small, KL renormalization 
cannot overshoot bare repulsion, and the bound state does not appear. The situation 
may change when we include momentum dependence of the interaction and non-circular 
nature of the pockets. In this last case, there appears infinite number of Aig harmonics, 
which all couple to each other, and in some cases one or several eigenfunctions may 
end up being attractive®^’®®. Besides, angle dependence generates d—wave and gf—wave 
harmonics, and some of eigenfunctions in these channels may also become attractive and 
compete with s—wave®^’®®. Still, however, in distinction to the isotropic case, there is 
no guarantee that “some" eigenfunction from either Aig, or Big, or i? 2 g, or A 2 g, will be 
attractive. In other words, a lattice system may well remain in the normal state down to 
T = 0. 


RPA-type approach, spin-mediated interaction 


How can we still get superconductivity in this situation? One way to proceed is to 
apply another ladder summation scheme - this time to series of renormalizations which 
transform a bare interaction into an irreducible particle-particle vertex. The leading 
terms in the series are KL terms, but full ladder series include inhnite set of higher- 
order terms. This computational procedure is often called random-phase approximation 
(RPA) by analogy with the analogous summation scheme to get a screened Coulomb 
interaction. I skip the details of the calculations (they can be found in, e.g., ^ 2,100 
and formally require npft(O) 3> Ilph{2kF) and Yiph{Q) 3> Aph{Q + ^kp)) and present 
the result: ladder summation gives rise to an irreducible pairing vertex in the form 
-/c;p, -p) = Tc{k - p)6a-ydi3s + Ts{k - p)Raa ■ ^135, where for k and p on the 
same pocket 


rc(o) = 


Ga 


r.(o) = 


Ga 


2 1 -\- G4npft(0) ^ 2 1 — G4np/i(0) 

and for k and p at different pockets, when k — p ^ Q 


(14) 


rc(Q) = 




2 1 -\- GzAphiQ) 


r.(g) = 




2 1 — GsIlph{Q) 


(15) 
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Figure 12: The plot of n(g) for a 2-pocket model with q along the zone diagonal. When q < 2kF, 
n(q) saturates, as it is expected for a 2D system with a circular Fermi surface. Note the 2kF 
cusp-like feature, which is the one-sided 2kF non-analyticity of n(g') in 2D. At larger q, n(g) 
gets larger and almost diverges at g ~ Q due to near-nesting. The inset shows the FS topology 
for which n(q) has been calculated. The arcs at the corners are parts of the electron pocket and 
the one in the center is the hole pocket. From [32,101]. 
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Re-expressing —k]p, —p) in terms of singlet and triplet components as 

^%,^5ik,-t,p, -p) = 

Tg—o(/c p) “ 1 “ 

rs=l(^ p) “t“ ^q:(5^/37) • 

(16) 


we obtain 


rs=o — - (Fc — SFs 
Fs=i = - (Fc -f F^) 


r _ in') = -— ( ^ ^ ^ 

’ 4 + G4np70) ^ 1 - G,Uph{0)) 

r _ roi = ^ ^ _ ^ ^ 

^ 4 '^i-G'4np7o) i + G4np,(o)y 

P _ ^ 1 , 3 \ 

4 Vl + GsRp^g) l-GsUph{Q)) 

F3=i(g) = — (r^^sivM ~ iTGsIvm) 

Let’s compare this result with what we obtained in the KL formalism. Focus on the 
singlet channel and expand in (18) to second order in (^ 3 ^ 4 . We have 

F,=o(0) ^ -7- (^ + 

~ —G 4 (1 -|- 0.5G4Flp/i(0)) 

F,=o(g) ~ - 7 - (^ + l-GsRp^g)) 

~ —G 3 (1 -|- 0.5G3Flp/i(g)) (19) 


Apart from the factor of 1/2 (which is the consequence of an approximate RPA scheme) 
Fs=o(0) is the same as irreducible vertex F®^, which we obtained in KL calculation in the 
previous section, and Fs=o(g) the same as T 12 By itself, this is not surprising, as in Fs=o 
we included the same particle-hole renormalization of the bare pairing interaction as in 
the KL formalism. 


I now look more closely at the spin-singlet components 


Fs=o(0) 

Fs=o(g) 


1 / G4 3^4 \ 

4^1-1- G4Flp/j(0) 1 — G4Flpft(0 )) 
1 ( Gs 3 G 3 

■4 ^ 1 - G^Ii^h{Q) 


(20) 
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For repulsive interaction, the charge contribution gets smaller when we add higher terms 
in G whereas spin contribution gets larger. A conventional recipe in this situation is to 
neglect all renormalizations in the charge channel and approximate rs=o with the sum 
of a constant and the interaction in the spin channel. The irreducible interaction in the 
s — channel is then 


rs=o(o) — rg=o(Q) — 

-^4 + Gs 3 / G 4 
4 4 l,l-G4np;,(0) 


1 


G3 \ 

GsUpi^iQ)) 


( 21 ) 


Like I said before, if G 4 np/i( 0 ) and GsUphiQ) are both small, G 4 — G 3 term is the 
largest and the pairing interaction is repulsive for G 4 > G 3 . However, we see that there 
is a way to overcome the initial repulsion: if GsUphiQ) > G 4 np/i( 0 ), one can imagine a 
situation when G^UphiQ) ~ 1 , and the correction term in ( 21 ) becomes large and positive 
and can overcome the negative hrst-order term. 

What does it mean from physics perspective? We found earlier that the condition 
GsUphiQ) = 1 signals an instability of a metal towards a SDW order with momentum Q. 
We don’t need the order to develop, but we need SDW fluctuations to be strong and to 
mediate pairing interaction between fermions. Once spin-mediated interaction exceeds 
bare repulsion, the irreducible pairing interaction in the corresponding channel becomes 
attractive. Notice in this regard that we need magnetic fluctuations to be peaked at large 
momentum transfer Q. If they are peaked at small momenta, np/j(0) exceeds 'n.ph{Q)-, 
and the interaction in the singlet channel remains repulsive. 


Spin-fluctuation approach 


What I just described is the main idea of the spin-fluctuation-mechanism of super¬ 
conductivity. The effective pairing interaction can be obtained either within 
or, using one of several advanced numerical methods developed over the last decade, or 
just introduced semi-phenomenologically. The semi-phenomenological model is called the 
spin-fermion modeT®^. Qnite often, interaction mediated by spin flnctuations also crit¬ 
ically affects single-fermion propagator (the Green’s fnnction), and this renormalization 
has to be included into the pairing problem. As another complication, the interaction 
mediated by soft spin flnctnations has a strong dynamical part dne to Landau damping 
- the decay of a spin fluctnation into a particle-hole pair. This dynamics also has to 
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be included into consideration, which makes the solution of the pairing problem near a 
magnetic instability quite involved theoretical problem. 

There are two crucial aspects of the spin-fluctuation approach^®^’^®^. First, magnetic 
fluctuations have to develop at energies much larger than the ones relevant for the pairing, 
typically at energies comparable to the bandwidth W. It is crucial for spin-fluctuation 
approach that SDW magnetism is the only instability which develops at such high en¬ 
ergies. There may be other instabilities (e.g., charge order), but the assumption is that 
they develop at small enough energies and can be captured within the low-energy model 
with spin fluctuations already presentSecond, spin-fluctuation approach is fun¬ 
damentally not a weak coupling approach. In the absence of nesting, Yiph{Q) and np/i(0) 
are generally of order l/hF, and 'n.ph{Q) is only larger numerically. Then the interaction 
G 3 must be of order W in order to get a strong magnetically-mediated component of the 
pairing interaction. 

One way to proceed in this situation is to introduce the spin-fermion model with 
static magnetic fluctuations built into it, and then assume that within this model the 
interaction between low-energy fermions g is smaller than W and do controlled low-energy 
analysis treating g/W as a small parameter^°^’^°^’^®^. There are several ways to make 
the assumptions g W and G ~ hF consistent with each other, e.g., if microscopic 
interaction has length Tq and T^kp/h 3> 1, then g is small in compared to G 

(Refs.^^^’^'^®). At the same time, the properties of the spin-fermion model do not seem 
to crucially depend on g/W ratio, so the hope is that, even if the actual g is of order 
W, the analysis based on expansion in g/W captures the essential physics of the pairing 
system behavior near a SDW instability in a metal. 


IV. INTERPLAY BETWEEN SDW MAGNETISM AND SUPERCONDUCTIV¬ 
ITY, PARQUET RG APPROACH 

I now return to weak coupling, where I have control over calculations, and ask the 
question whether one can still get an attraction in at least one pairing channel despite 
that G 4 > G 3 , i.e., the bare pairing interaction is repulsive in all channels. The answer 
is, actually, yes, it is possible, but under a special condition that Yiph{Q) is singular and 
diverges logarithmically at zero frequency or zero temperature, in the same way as the 
particle-particle bubble npp(O). This condition is satisfied exactly when there is a perfect 



nesting between fermionic excitations separated by Q. For Fe-pnictides, it implies that 
hole and electron FSs perfectly match each other when one is shifted by Q. 

I show below that 'n.ph{Q) and npp(O) do have exactly the same logarithmic singularity 
at perfect nesting. At the moment, let’s take this for granted and compare the relevant 
scales. First, no fluctuations develop at energies/temperatures of order W because at 
such high scales the logarithmical behavior of Ylpp and liph is not yet developed and 
both bubbles scale as 1/W. At weak coupling G/W « 1, hence corrections to bare 
vertices are small at these energies. Second, we know that the pairing vertex evolves at 
(G 3 —G 4 )npp( 0 ) ~ 0 ( 1 ), and that corrections to the bare irreducible pairing vertex become 
of order one when Gz^phiQ) ~ 0(1). But we also know from, e.g., (15) that at the same 
scale the SDW vertex begins to evolve. Moreover other inter-pocket interactions, which we 
didn’t include so far: density-density and exchange interactions (which here and below 
we label as Gi and G 2 , respectively) also start evolving because their renormalization 
involves terms GiIlph{Q) and G 2 llph{Q), which also become of 0 ( 1 ), provided that all 
bare interactions are of the same order. Once Gi^ 2 ^ph{Q) becomes of order one, the 
renormalization of G 3 by Oi and G 2 interactions also becomes relevant. The bottom 
line here is that renormalization of all interactions become relevant at the same scale 
where GiUphiQ) ~ Gjnpp(O) ~ 1 . At this scale we can expect superconductivity, if the 
corrections to G 4 — G 3 overcome the sign of the pairing interaction, and we also we can 
expect an instability towards SDW and, possibly, towards some other order. The issue 
then is whether it is possible to construct a rigorous description of the system behavior in 
the situation when all couplings are small compared to W, but GiIlph{Q) and Ginpp(O) 
are of order one. The answer is yes, and the corresponding procedure is called a parquet 
renormalization group (pRG). 

The pRG is a controlled weak coupling approach. It assumes that no correlations 
develop at energies comparable to the bandwidth, but that there are several competing 
orders whose fluctuations develop simultaneously at smaller energies. Superconductivity 
is one of them, others include SDW and potential charge-density-wave (GDW), nematic 
and other orders. The pRG approach treats superconductivity, SDW, GDW and other 
potential instabilities on equal footings. Gorrelations in each channel grow up with similar 
speed, and fluctuations in one channel affect the fluctuations in the other channel and 
vise versa. For superconductivity, once the corrections to the pairing vertex become of 
order one, and there is a potential to convert initial repulsion into an attraction. We 
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know that second-order contribution to the pairing vertex from SDW channel works in 
the right direction, and one may expect that higher-order corrections continue pushing 
the pairing interaction towards an attraction. However even if attraction develops, there 
is no guarantee that the system will actually undergo a SC transition because it is entire 
possible that SDW instability comes before SC instability. 

The pRG approach addresses both of these issues. It can be also applied to a more 
realistic case of non-perfect nesting if deviations from nesting are small in the sense that 
there exists a wide range of energies where liphiQ) and npp(O) are approximately equal. 
Below some energy scale, cuo, the logarithmical singularity in Ylph{Q) is cut. If this scale 
is smaller than the one at which the leading instability occurs, a deviation from a perfect 
nesting is an irrelevant perturbation. If it is larger, then pRG runs up to cuo, and at 
smaller energies only SC channel continues to evolve in BCS fashion. 

There also exists a well-developed numerical computational procedure called functional 
RG Its advantage is that it is not restricted to a small number of patches 

and captures the evolution of the interactions in various channels even if the interactions 
depend on the angles along the FS. The “price" one has to pay is the reduction in the 
control over calculations - fRG includes both leading and subleading logarithmical terms. 
If only logarithmical terms are left, the angle dependencies of the interactions do not 
evolve in the process of RG flow, only the overall magnitude changes^®® So far, the results 
of fRG and pRG analysis for various systems fully agree. Below I focus on the pRG 
approach. For the thorough tutorial on the RG technique, see Ref.^'^®. In the discussion 
below and in Sec. 8.5 I follow Refs. [32,101]. 

A. Parquet Renormalization Group: The Basics 

I recall that in Fe-pnictides a bubble with momentum transfer Q contains one hole 
(c) and one electron (f) propagator, and at perfect nesting the dispersions of holes and 
electrons are just opposite, ec{k) = —ef{k + Q). The particle-hole and particle-particle 
bubbles are 


n,,(0) = 

r d'^k dij .1 

n„.( 0 ) = » J u:)Gk{Q + k,uj), (22) 

where 
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= ^‘^•f+iSsgn{uj) ' Substituting into Eq. 22 and using edk) = —Sf{k + Q) one can 
easily make sure that the two expressions in Eq. 22 are identical. Evaluating the integrals 
we obtain 

npp(O) = Tlph{Q) = NqL + ... (23) 

where Nq = m/2nh‘^ is the 2D density of states, 



E is a typical energy of an external fermion, and the dots stand for non-logarithmic terms. 
The factor 1/2 is specihc to the pocket model and accounts for the fact that for small 
pocket sizes, the logarithm comes from integration over positive energies W > E > Ep. 
At non-perfect nesting, the particle-particle channel is still logarithmic, but the particle- 
hole channel gets cut by the energy difference {SE) associated with the nesting mismatch, 
such that 


n,»W) = iV„log-^D^ (25) 

The main idea of pRG (as of any RG procedure) is to consider E as a running variable, 
assume that initial E is comparable to W and Gi log = GiL is small, calculate the 
renormalizations of all couplings by fermions with energies larger than E, and hnd how 
the couplings evolve as E approaches the region where G^L = 0(1). 

This procedure can be carried out already in BGS theory, because Gooper renormal¬ 
izations are logarithmical. For an isotropic system, the evolution of the interaction Ui in 
a channel with angular momentum I due to Gooper renormalization can be expressed in 
RG treatment as an equation for the running coupling Ui{L) 

'EM ^ .N„{U,(L)f. (26) 


The solution of (26) is 


Ui{L) = 


Ui 


(27) 


1 + UiN„L 

Similar formulas can be obtained in lattice systems when there are no competing instabil¬ 
ities, i.e., only renormalizations in the pairing channel are relevant. For example, in the 
two-pocket model for the pnictides, the equations for the vertices Vhh{,L) = —Gi{L) and 
Vhe{,L) = —G^i^L), Eqs. (8), can be reproduced by solving the two coupled RG equations 

dGs{L) 


dL 

dG^{L) 

dL 


= -2NoGsiL)G^iL) 

= -/V„ + (Gt(L) f) 


(28) 
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with boundary conditions Gi{L = 0) = G 4 , G 3 (L = 0) = G 3 . The set can be factorized 
by introducing Ga{L) = G^{L) + Gi{L) and Gb{L) = G 4 {L) — G 2 ,{L) to 


= -/V„ (0^(1))", = -N„ (Gfl(L))^ 

dL dL 


(29) 


The solution of the set yields 

Ga{L) = G4{L) + G^{L) = 
Gb{L) = GM - G,{L) = 


G3 + G4 

1 + iVoT(G3 + G4) 
G4 — G3 


(30) 


1 + NqL{G 4 — G 3 ) 

Solving this set and using Thh{L) = —G 4 {L), Tfie{L) = —G^i^L), we reproduce ( 8 ). This 
returns us to the same issue as we had before, namely if G 4 > G 3 , the fully renormalized 
pairing interaction does not diverge at any L and in fact decays as L increases: G 4 {L) 
decays as 1 /L and G^i^L) decays even faster, as 1 /L^. 

I now consider how things change when 'n.ph{Q) is also logarithmical and the renormal¬ 
izations in the particle-hole channel have to be included on equal footings with renormal¬ 
izations in the particle-particle channel. 


B. pRG in a 2-pocket model 

Because two types of renormalizations are relevant, we need to include into consid¬ 
eration all vertices with either small total momentum or with momentum transfer near 
Q i.e., use the full low-energy Hamiltonian of Eq. (2). There are couplings G3 and G4 
which are directly relevant for superconductivity, and also the couplings Gi and G 2 for 
density-density and exchange interaction between hole and electron pockets, respectively. 
These are shown in Fig. 8 . 

The strategy to obtain one-loop pRG equations, suitable to our case, is the following: 
One has to start with perturbation theory and obtain the variation of each full vertex 
5Gi to order GiGjL. Then one has to replace 6Gi/L by dGi{L)/dL and also replace 
GiGj in the r.h.s. by Gi{L)Gj{L). The result is the set of coupled differential equations 
for dGi{L)/dL whose right sides are given by bilinear combinations of Gi{L)Gj{L). The 
procedure may look a bit formal, but one can rigorously prove that it is equivalent to 
summing up series of corrections to G* in powers of GiL, neglecting corrections terms 
with higher powers of Gi than of L. One can go further and collecting correction terms 
of order GiGjGkL. This is called 2-loop order, and 2-loop terms give contributions of 
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order (G(L))^ to the right side of the equations for dGi{L)/dL. 2-loop calculations are, 
however, quite involved^^^ and have not been re-checked. Below I only consider 1-loop 
pRG equations. 



ir • .ix 


93^ 





Figure 13: The pRG diagrams to one loop order, which contribute to the parquet flow of gi, g 2 , 
gs and g^ vertices. From [32,101]. 

The corrections to all four couplings are shown in Fig.13. Evaluating the integrals 
and following the recipe we obtain 

91 = 91 + 93 

92 = 25'2(5'i ~ 92 ) 

93 = ‘^93{‘^9i — 92 — 9i) 

94: = —93 ~ 94 

(31) 

where we introduced gi = gi{L) = Gi{L)NQ and cji = dgi/dL 

We note that the renormalizations of g 4 are still only in the Cooper channel and causes 
g 4 to reduce. But for g^ we now have a counter-term from gi, which pushes g 3 up. And 
the gi term is in turn pushed up by 5 ^ 3 . Thus already at this stage one can qualitatively 
expect g^ to eventually get larger. Fig 14 shows the solution of (31)- the flow of the four 
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couplings for this model. We see that, even if is initially smaller than g 4 , it flows up 
with increasing L, while 5^4 flows to smaller values. At some L = Lq, g^ crosses g^, and 
at larger L the pairing interaction g 4 — g^ becomes negative (i.e., attractive). In other 
words, in the process of pRG flow, the system self-generates attractive pairing interaction. 
I remind that the attraction appears in the channel. The pairing interaction in s’*"’*' 
channel: g^ + g 4 remains positive (repulsive) despite that g 4 eventually changes sign and 
becomes negative. It is essential that for L ^ Lq the renormalized gi are still of the same 
order as bare couplings, i.e., are still small, and the calculations are fully under control. 
In other words, the sign change of the pairing interaction is a solid result, and higher-loop 
corrections may only slightly shift the value of Lq when it happens. 

At some larger L = Lc, the couplings diverge, signaling the instability towards an 
ordered state (which one I discuss later). One-loop pRG is valid "almost" all the way 
to the instability, up to Lc ~ ~ when the renormalized gi become of order one. 

At smaller distances from Lc higher-loop corrections become relevant. It is very unlikely, 
however, that these corrections will change the physics in any significant way. 

The sign change of the pairing interaction can be detected also if the nesting is not 
perfect and 'n.ph{Q) does not behave exactly in the same way as npp(O). The full treatment 
of this case is quite involved. For illustrative purposes I follow the approach first proposed 
in Ref.^^^ and measure the non-equivalence between npp(O) and liphiQ) by introducing 
a phenomenological parameter di = Yiph{Q)/^ppiS^) and treat di as an L— independent 
constant 0 < di < 1, independent on L. This is indeed an approximation, but it is at least 
partly justified by our earlier observation that the most relevant effect for the pairing is 
the sign change of g 4 — g^ at some scale Lq, and around this scale di is not expected to 
have strong dependence on L. The case di = 1 corresponds to perfect nesting, and the 
case di = 0 implies that particle-hole channel is irrelevant, in which case, I remind, (74 —(73 
remains positive for all L. 

The pRG equations for arbitrary di are straightforwardly obtained using the same 
strategy as in the derivation of (31), and the result 


9i = di{9i + 9z) 
g2 = 2di5f2(fi'i ~ 92 ) 

93 = ‘^dig3{2gi — g2) — 25 ( 35(4 

94 = - 9 I - 9 I 

(32) 
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Figure 14: The flow of dimensionless couplings 51,2,3,4- 53 grows and eventually crosses 54 , which 
becomes negative at a large enough RG scale. From [112]. 
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RG scale L 
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RG scale L 


Figure 15: The flow of ratio of couplings gzldi and g^lgi for different nesting parameters 
di = 1(a), di = 0.3(b),di = 0.05(c). All cases are qualitatively similar in that g^/gi eventually 
crosses gijgi. The smaller is the nesting parameter, the ‘later’ is this crossing. If di = 0, this 
crossing will never happen and 54 > 53 for all L. 


In Fig 15 I show the behavior of the couplings for representative 0 < di < 1. Like 
before, I take bare value of 5^4 to be larger than the bare i.e., at high energies the 
pairing interaction is repulsive. This hgure and analytical consideration shows that for 
any non-zero di the behavior is qualitatively the same as for perfect nesting, i.e., at some 
Lo < Lc the running couplings g^, and 5^4 cross, and for larger L (smaller energies) pairing 
interaction in s'*"” channel becomes attractive. The only effect of making di smaller is the 
increase in the value of Lq. Still, for sufficiently small bare couplings, the range where the 
pairing interaction changes sign is fully under control in one-loop pRG theory. 

A way to see analytically that g^ — g^ changes sign and becomes positive is to consider 
the system behavior near L = L^. and make sure that in this region gz> g a- One can easily 


make sure that all couplings diverge at Lc, and their ratios tend to some constant values 
(see discussion around Eq. (42) below for more detail). Introducing g 2 = agi^g^ = bgi, 


and gA = cgi, and substituting into (32) we hnd an algebraic set of equations for a, b, and c. 


Solving the set, we hnd that b = 


-and c = ^(3 — 6 ^). The negative sign 


of c and positive sign of b, combined with the fact that gi dehnitely increases under the 


how and surely remains positive, imply that near Lc, g 4 is negative, while g^ is positive 
(this is also evident from the Fig 15). Obviously then, g^ and gA must cross at some 
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Lq < Lc- 


The reason for the sign change of the pairing interaction is clear from the structure of 
the pRG equation for the r.h.s. of which contains the term Adig^g^, which pushes g^ up. 
We know from second-order KL calculation that the upward renormalization of g^ comes 
from the magnetic channel and can be roughly viewed as the contribution from spin- 
mediated part of effective fermion-fermion interaction. Not surprisingly, we will see below 
that gi does, indeed, contribute to the SDW vertex. From this perspective, the physics 
of the attraction in pRG (or in fRG, which brings in the same conclusions as pRG) and 
in spin-fermion model is the same: magnetic fluctuations push inter-pocket/inter-patch 
interaction up, and below some energy scale the renormalized inter-pocket/inter-patch 
interaction becomes larger than repulsive intra-pocket/intra-patch interaction. 

There is, however, one important difference between the RG description and the de¬ 
scription in terms of spin-fermion model. In the spin-fermion model, magnetic fluctuations 
are strong, but the system is assumed to be at some distance away from an SDW insta¬ 
bility. In this situation, SG instability dehnitely comes ahead of SDW magnetism. There 
may be other instabilities produced by strong spin fluctuations, like 
which compete with SG and, by construction, also occur before SDW order sets in. 

In RG treatment (pRG or fRG), SDW magnetism and SG instability (and other po¬ 
tential instabilities) compete with each other, and which one develops hrst needs to be 
analyzed. So far, we only found that SG vertex changes sign and becomes attractive. 
But we do not know whether superconductivity is the leading instability, or some other 
instability comes hrst. This is what we will study next. The key issue, indeed, is whether 
superconductivity can come ahead of SDW magnetism, whose huctuations helped convert 
repulsion in the pairing channel into an attraction. 


V. COMPETITION BETWEEN DENSITY WAVE ORDERS AND SUPERCON¬ 
DUCTIVITY 

Thus far, we identihed an instability in a particular channel with the appearance of a 
pole in the upper frequency half-plane in the corresponding vertex - the vertex with zero 
total momentum in the case of SC instability, and the vertex with the total momentum 
Q in the case of SDW instability. Since our goal is to address the competition between 
these states, it is actually advantageous to use a slightly different approach: introduce 
all potentially relevant fluctuating helds, use them to decouple 4-fermion terms into a set 
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of terms containing two fermions and a fluctuating field, compute the renormalization of 
these “three-legged" vertices and use these renormalized vertices to obtain the suscepti¬ 
bilities in various channels and check which one is the strongest. We will see that the 
renormalized vertices in different channels (most notably, SDW and SC) do diverge near 
Lc, but with different exponents. The leading instability will be in the channel for which 
the exponent is the largest. There is one caveat in this approach — for a divergence of 
the susceptibility the exponent for the vertex should be larger than 1/2 (Ref.®'^), but we 
will see below that this condition is satisfied, at least for the leading instability. 

A. Two pocket model 

Let us see how it works for a two-pocket model. There are two particle-particle three 
legged vertices Vh,e as shown in Fig 16. To obtain the flow of these vertices, i.e., rf^(L) 
I assume that external fermions and a fluctuating field have energies comparable to some 
E (i.e.,L = log A/E) and collect contributions from all fermions with energies larger than 
E. To do this with logarithmical accuracy I write all possible diagrams, choose a particle- 
particle cross-section at the smallest internal energy E' > E and sum up all contributions 
to the left and to the right of this cross-section, as shown in Fig 17. The sum of all 
contributions to the left of the cross-section gives the three legged vertex at energy E' (or 
L' = log A/E'), and the sum of all contributions to the right of the cross-section gives the 
interaction gi at energy L'. The integration over the remaining cross-section gives dL' 
(with our normalization of Qi), and the equation for, e.g., Th{L) becomes 

rf (L) = f dV (rf^(L')j4(r') + rf^(L')j3(r')) (33) 

Differentiating over the upper limit, we obtain differential equation for dTf^^{L)/dL whose 
r.h.s. contains Ff^(L) and gz,i{,L) at the same scale L. 

Collecting the contributions for Ff‘"(L) an T^^[L) we obtain 

^ h 9i + ^ e 93 

rf (74 + Ti^93 

(34) 


dir 

dL 

d^ 

dL 


or 
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pSDW 


pCDW 


Figure 16: Superconducting and density-wave three-leg vertices. Divergence of any of these 

vertices indicates that the system is likely to be unstable to the corresponding order. F^^ are 

superconducting vertices, is SDW vertex and is CDW vertex. From [112]. 
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Figure 17: The flow diagrams for the effective vertices: SC vertex (top two) and SDW vertex 
(bottom). The couplings gi's here are running couplings in RG sense. From [32,101]. 
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(35) 


d^++ 

dL 

dT+_ 

dL 


(fi'4 + fi'3)r++ 
(^4 -^3)r+- 


where r++ = Ff'" + Ff*" and F+_ = Ff'" — Ff*". The first one is for s’*""'' pairing, the 
second is for s’*"” pairing. We have seen in the previous section that the running couplings 
5 f 3 4 diverge at some critical RG scale Lc- The flow equation near is in the form g ^ 
hence 


9i = 



(36) 


Substituting this into Eq. 35 and solving the differential equation for F we hnd that the 
two SC three legged vertices behave as 


Fs++ oc 




1 


F.+- 


1 

{Lc - L)“3-«4’ 


(37) 


The requirement for the divergence of F 5 +- is 0:3 > 0:4, which is obviously the same as 
gs > fi '4 (see (36)). 

I follow the same procedure for an SDW vertex F^^'^ = T^^^dai 3 - I introduce a 
particle-hole vertex with momentum transfer Q and spin factor cTq,^, as shown in Fig 16, 
and obtain the equation for dT^^^{L)/dL in the same way as we did for SC vertices. We 
obtain (see Fig. 17) 


^YSdw 

dL 


diigi + g3)T^^^ 


(38) 


Using Eq. 36 and following the same steps as above we obtain at L ~ Lc 


^SDW 


(Lc - L)'^lCl+a3) 

For CDW vertex (the one with the overall factor 6 ai 3 instead dag), the flow equation is 


^^CDW 

dL 


d,{g, + g,-2gs-2g2)T^'^^ 

di{gi - gs- 2g2)T^^^ 


(40) 


Using the same procedure as before we obtain 

YCDw ^ ^ _ 

— JJ^di{ai-a3-2a2) 


(41) 
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The exponents a* can be easily found by plugging in the asymptotic forms in Eq. 36 
into the RG equations. This gives the following set of non linear algebraic equations in ctj 

Ui = T ctg) 

02 = 2,diO!2{o!i — 02 ) 

03 = 2 , diC ( s ( 2 , o'i — 02) — 20304 

04 = —O3 — O4 

(42) 

Consider first the case of perfect nesting, di = 1. The solution of the set of equations 
is Oi = |, 02 = 0, 03 = ^ and 04 = — Combining o’s, we hnd that the exponents for 
superconducting and spin density wave instabilities and positive and equal: 


Ois± = O 3 — O 4 — 


OiSDW = Oi + O 3 — 


i + Vs 

6 

i + Vs 

6 


^ 0.539 
^ 0.539 


(43) 


while the exponent for CDW and s + + vertices are negative 

1 - ys 


OiCDW — Oi + O 3 — 


Os+_|_ — —O 3 — O 4 — 


6 

1 - 75 

"6 


0.206 

-0.206 


(44) 


We see that the superconducting (s’*"”) and SDW channels have equal susceptibilities 
in this approximation, while CDW channel is not a competitor. 

The analysis can be extended to di < 1. I dehne fd = q; 4 /q;i, 7 = a^jo-x and obtain 


7^ = 


^164 - 4df + 4 + 2 - d? 


df 


1(3-74 


2 

1 1 


0.1 — 


(45) 


di 1 -4 7 ^ 

In FiglS I plot as± = — 0 : 4 , asnw = cti + 0 : 3 , and acDW = cni — 0 : 3 , We clearly see 

that (i) CDW channel is never a competitor, and (ii) as di decreases (the nesting gets 
worse), the pairing vertex diverges with a higher exponent that SDW channel, hence s’*"” 
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superconductivity becomes the leading instability, overshooting the channel which helped 
SC vertex to change sign in the hrst place. 

In real systems, pRG equations are only valid up to some distance from the instability 
at Lc- Very near Lc three-dimensional effects, corrections from higher-loop orders and 
other perturbations likely affect the flow of the couplings. Besides, in pocket models, 
the pRG eqnations are only valid for E between the bandwidth W and the Fermi energy 
Ep. At E < Ep, internal momenta in the diagrams, which acconnt for the flow of the 
conplings, become smaller than external kp, and the renormalization of gi start depending 
on the interplay between all fonr external momenta in the vertices^^’^^®. The calcnlation 
of the flow in this case is technically more involved, but the result is physically transparent 
- SDW and s’*"” SC channels stop talking to each other, and the vertex evolves according 
to Eqs. (37) and (38), with gi taken at the scale Ep (or Lp = \ogA/Ep). If Lp > 
the presence of the scale set by the Fermi energy is irrelevant, bnt if Lp < Lc (which is 
the case for the Fe-pnictides because superconducting Tc and magnetic Tsdw are much 
smaller than Ep), then one should stop pRG flow at Lp^. At perfect nesting, the SDW 
combination gi + g^ is larger than s’*"” combination g^ — g 4 at any L < Lc, hence SDW 
channel wins, and the leading instability npon cooling down the system is towards a 
SDW order. At non-zero doping, flph{Q) is cnt by a deviation from nesting, what in onr 
language implies that di < 1. If bare g^ and g^ are not to far apart, there exists a critical 
di at which g^ — g^ crosses di{gi + gz) at Lp, and at larger di the crossing occurs before 
Lp. In this sitnation, s^~ SG becomes the leading instability npon cooling off the system. 

The comparison between different channels can be further extended by considering 
cnrrent SDW and GDW vertices (imaginary and F'"^'^) and so on. I will not dwell 

into this issue. 

Before moving on, I need to clarify one more point. So far we fonnd that the vertices 
F‘^‘" and diverge and compared the exponents. However, to actnally analyze the 

instability in a particular channel one has to compnte flnctnation correction to suscepti¬ 
bility 


x)AL) ~ / d^k{rf If, oc fdE {r{L')f (46) 

where Hj is either H"^^^ = Tfph or H'^'" = Tfpp (see Fig 19) 

The fully renormalized susceptibility in a given channel is 

{x'(L)A = ri-x)m (47) 
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Figure 18: Exponents (as±, asDW and acow) for different values of the nesting parameter di 
calculated near the critical RG scale, where the couplings diverge. The state with the largest ex¬ 
ponent wins. SDW and SC are degenerate when di = 1 (perfect nesting) and superconductivity 
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wins for all other values of di. CDW is not a competitor for all values of di. 









Figure 19: (Left) The fluctuation correction to SC pairing susceptibility. (Right) The fluctuation 
correction to SDW susceptibility. 
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where Tq is some bare value of order one. The true instability occurs at L* when = 

Tq. At weak coupling, the critical L* is close to Lc, and, indeed, the instability occurs 
hrst in the channel with the largest exponent for T*. However, we need x}i{L) to diverge 
at Lc, otherwise there will no instability at weak coupling®’^. This requirement sets the 
condition that the exponent for the corresponding T must be larger than 1/2. Fortunately, 
this condition is satished in the two-pocket model. For di = 1, this is evident from (43). 
For di < 1, the exponent for the SC channel only increases, while the one in SDW channel 
decreases but still remains larger than 1/2 as it is evidenced from FiglS where I plotted 
the exponents for SC and SDW vertices as a function of di. In the limit di —)• 0, 


1 di 
asDW ~ 2 ^ T 


(48) 


The fact that both asc and asow are larger than 1/2 implies that in Landau-Ginzburg 
expansion in powers of SC and SDW order parameters (A and M, respectively), not only 
the prefactor for changes sign at T^, but also the prefactor for term changes sign 
and becomes negative below some < T^. This brings in the possibility that at low T 
SC and SDW orders co-exist. The issue of the co-existence, however, requires a careful 
analysis of the interplay of prefactors for fourth order terms A^, and M^A^. I do not 
discuss this specihc issue. For details see^^®’^^®. 


1. Multi-pocket models 

The interplay between SDW and SC vertices is more involved in more realistic multi¬ 
pocket models Fe-pnictides, with several electron and hole pockets. I recall that weakly 
doped Fe-pnictides have 2 electron pockets and 2-3 hole pockets. In multi-pocket models 
one needs to introduce a larger number of intra-and inter-pocket interactions and analyze 
the flow of all couplings to decide which instability is the leading one. This does not 
provide any new physics compared to what we have discussed, but in several cases the 
interplay between SC and SDW instabilities becomes such that superconductivity wins 
already at perfect nesting. In particular, in 3-pocket models (two electron pockets and one 
hole pockets) the exponent for the SC vertex gets larger than the exponent for the SDW 
vertex already at di = 1. I show the flow of SC an SDW couplings for 3-pocket model 
in Fig.20. Once di becomes smaller than one, SC channel wins even bigger compared to 
SDW channel. 
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Figure 20: The flow of the SC and SDW vertices with the RG scale. Both diverge at a critical 
scale, Lc, but the SC vertex diverges stronger. From Ref. [108]. 
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Superconductivity right at zero doping has been detected in several Fe-pnictides, like 
LaOFeAs and LiFeAs, and it is quite possible that this is at least partly due to the specihcs 
of pRG flow. 

B. Summary of the pRG approach 

I now summarize the key points of the pRG approach 

• The SG vertex starts out as repulsive, but it eventually changes sign at some RG 
scale (To). This happens due to the "push" from SDW channel, which rives rise to 
upward renormalization of the inter-pocket interaction g^. 

• Both SDW and SG vertices diverge at RG scale Lc which is larger than Lq. The 
leading instability is in the channel whose vertex diverges with a larger exponent. 
At perfect nesting, SDW instability occurs first in 2-pocket model, however in some 
multi-pocket models SG vertex has a larger exponent that the SDW vertex and SG 
becomes the leading instability. 

• Deviations from perfect nesting (quantihed by di < 1) act against SDW order by 
reducing the corresponding exponent. At snfficiently small di SG instability becomes 
the leading one. 

• The necessary condition for the instability is the diverges of the flnctnating compo¬ 
nent of the susceptibility. This sets up a condition a > ll2, where a is the exponent 
for the corresponding vertex. For the leading instability, we found a > 1/2 in all 
cases. For the snbleading instability, a can be either larger or smaller than 1/2. 
This affects potential co-existence of the leading and snbleading orders at a lower 
T. 

VI. SDW MAGNETISM AND NEMATIC ORDER 

For this section, I assnme that we are in the range of parameters/dopings, where SDW 
instability comes hrst, and consider (i) what kind of SDW order emerges and (ii) the 
interplay between breaking of 0(3) spin-rotational symmetry and breaking of a discrete 
O 4 symmetry of rotations on a tetragonal lattice. I consider these two issues one after the 
other. In the discussions in this section I follow Refs. [63,64,120]. 
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A. Selection of SDW order 


I return to the model I started with, but now with interactions renormalized by pRG 
contributions from energies lager than Ep. The only necessary extension we need to 
make is we need to consider two electron pockets, one at (0,7r) and another at (7r,0) in 
the unfolded Brillouin zone (see Fig.21). To make presentation more simple, we consider 
only one hole pocket, centered at (0,0). The extension to two (or three) hole pockets is 
straightforward, but requites care and in some cases leads to new states^"^^’^^® 

We need to be a bit more precise and include the ellipticity of electron pockets. Ac¬ 
cordingly, we approximate dispersions of fermions near hole and electron pockets by 
£r,k = eo - ^ - fi = -e, £x,k+Qi = e - 5o + 82 cos 29, £y,k+Q 2 = e - 5o - 82 cos 29, 
where m* denotes the band masses, Sq is the offset energy, /i is the chemical potential, 
^0 = 2/i, 82 = Somirrix — niy) / {2mxmy) , and 9 = ky /■ 

I shift the momenta of the fermions near the X and Y Fermi pockets by Qi and Q 2 , 
respectively, i.e. £x,k+Qi -t £x,k, ^y,k+Q 2 -t Wk- 

This model has eight fermionic interactions (with the same structure as in a 2-pocket 
model, but now there are four different inter-and intra-pocket interactions involving the 
two electron pockets). These interactions can be decomposed into the spin density-wave 
(SDW), the charge density-wave (CDW) and the pairing channels. For magnetism, I keep 
only the interactions in the spin channel with momenta near Qi and Q 2 . This reduces 
the interacting Hamiltonian to 

'kiiat xl?spin ^ q (^ 9 ) 

i,q 

where Sj q = Y,k Cr,k+qa^«.sD,k /3 is the electronic spin operator, with Pauli matrices cTq,^. 
The coupling Ugpin is the combination of density-density and pair-hopping interactions 
between hole and electron states {gi and terms in the same notations as in previous 
two Sections). 


9lCp,aCr,aCx,i3Cx,0 = —^Cp^^(^al3Cx,i3 ' C^x jYT,S + (' ' ' ) 

fi'3Cr,«Cx,oCp,/3Cx,/3 = —^Cp,a^al3Cx,l3 ' Cx,'y^jYr,S + (' ' ' ) (50) 

where the dots stand for the terms with 5a,g8^^S) which only contribute to the CDW 
channel. Combining the two contributions for the SDW channel, I hnd gspm = gi + gs, 
as in ( 6 ). Once ggpin exceeds some critical value (which gets smaller when and 82 
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<A"x> = <Av> <A"x>><Av> 

<Ax> = <Av>=0 <A^ = <A,>=0 


Figure 21: (Color online.) (upper panel) The band-structure with a circular hole pocket at 
r and two electron pockets at X and Y. The Brillouin zone contains one Fe atom, (lower 
panels) Static magnetic susceptibility Xq across the Brillouin zone for different temperatures. 
At high temperatures, fluctuations near the two stripe magnetic ordering vectors are equally 
strong, (A^) = (Ay). Above the magnetic ordering temperature T/v but below the Ising- 
nematic ordering temperature T^, fluctuations associated with one of the stripe states become 
stronger (in the figure, (A^) > (Ay)) and the tetragonal symmetry is broken inside the unit 
cell. Stronger fluctuations around one ordering vector yield stronger intensity and narrower 
peaks. 
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decrease), static magnetic susceptibility diverges at (0,7r) and (tt, 0), and the system 
develops long-range magnetic order. An excitonic-type SDW instability in Fe-pnictides, 
resulting from the interaction between hole and electron pockets, has been considered by 
several authors^^®’^^^’^^^^^^®. 

My strategy is the following: I introduce the two bosonic helds A(x,y) oc 
EkC^. ka^o/ 3 C(x,y),k /3 for the collective magnetic degrees of freedom, use Hubbard- 
Stratonovich transformation to get rid of the terms in (49) with four fermions, integrate 
out the fermions, and obtain a Ginzburg-Landau (GL) action for and Ay. I then 
analyze this action in saddle-point approximation and show that one of the magnetic 
order parameters - either (Ay) or (Ay) - becomes non-zero in the magnetically ordered 
state. This leads to stripe-type SDW order in which spins are ordered ferromagnetically 
in one direction and antiferromagnetically in the other, i.e. the ordering momentum is 
either (tt, 0) or (0,7r). I then show that another state, in which (Ax) or (Ay) emerge 
simultaneously, may occur at a higher doping^^. The same tendency occurs in systems 
like Ba(Fei_a;Mna,) 2 As 2 , where the local Mn moments interact with the Fe conduction 
electrons^®. 

1. The action in terms of A.x and Ay 

A straightforward way to obtain the action in terms of (Ax) and (Ay) is to start with 
the fermionic Hamiltonian TL = TLq + Ttmt and write the partition function as the integral 
over Grassmann variables: 

Z oc y dci^\f^dcl y.e~^^ (51) 

and then decouple the quartic term in fermionic operators using the Hubbard-Stratonovich 
transformation: 



where, in our case, x = Sio = Hk<^vkoi'^al3^{x,Y),'kh V ~ ^(x,y)- One can then inte¬ 
grate Eq. (51) over fermionic variables using the fact that after the Hubbard-Stratonovich 
transformation the effective action becomes quadratic with respect to the fermionic op¬ 
erators. The result of the integration is recast back into the exponent and the partition 
function is expressed as: 
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Z oc J (53) 

If relevant Ay and Ay are small, which I assume to hold even if the magnetic tran¬ 
sition is first-order (I present the conditions on the parameters below), one can expand 
S'eflf [Ax, Ay] in powers of Ax and Ay and obtain the Ginzburg-Landau type of action 
for the order parameters Ax, Ay. For uniform Aj, the most generic form of S'eflf [ Ax, Ay] 

is 

S'eflf [Ax, Ay] = Tq (^A^ -|- Ay^ -I- - (^A^ -I- Ay^ 

+n(Ax-Ay)^ (54) 

Carrying out this procedure, one obtains the coefficients tq, u, g, and v in terms of the 
non-interacting fermionic propagators convoluted with Pauli matrices. The coefficient v 
vanishes in our model because of the anti-commutation property of the Pauli matrices: 

= 0 for i ^ j. To get a non-zero v, one needs to include direct interactions be¬ 

tween the two electron pockets^^°. The other three prefactors are expressed via fermionic 
propagators Gj^ = iun — ■Cj.k as 

=-^ ^ / ^r,fcGx,A: 

fi'spin 

W = - / Gp {Gx,k + Gy^kf' 

Z Jk ’ 

g = J^Gr,kiGx,k - Gy^kf (55) 

where = TY,nI and k = (k, a;„), with momentum k and Matsubara frequency 
Un = (2n -|- 1) ttT. Similar coefficients were found in Ref.^^°, which focused on the mag¬ 
netic instabilities in a two-band model. Near Tx,o one can expand vq as vq = a(T — Tx,o), 
with a > 0. Evaluating the integrals with the products of the Green’s functions, we obtain 

7C(3)iV^ 

7 / ~ - 

47r2T2 

g ^ 0.024m (56) 

for 6m T/sq 1. The crucial result for our consideration is that g is positive for any 
non-zero ellipticity. 

The action S'eflf is exact and includes all fluctuations of the two bosonic helds. Fluctu¬ 
ations need to be included for the analysis of a potential nematic order (see below), 
but the type of SDW can be analyzed already in the mean-field approximation (see 
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Rgfg63,i20 justification.) Solving for the minimum of S'efr[Ax, Ay] in Eq. (54), we 
find that, when (7 = 0 , the ground state has a huge degeneracy because any conhgura- 
tion A = (Ax)e*‘^i''' + (Ay) with (Ax)^ + (Ay)^ = —r^ju minimizes ^efr- A 

non-zero g gives rise to the additional coupling 2gAx^Y, which breaks this degeneracy. 
For a positive g, this term favors the states in which only one order parameter has a 
nonzero value, i.e. conhgurations with either (Ax) 7 ^ 0 or (Ax) 7 ^ 0 , but not both. 
These are stripe phases, in which spins order ferromagnetically along one direction and 
antiferromagnetically along the other one. 

For larger dopings, recent calculations^^ have shown that g may change sign and become 
negative. Then the SDW phase does not break C 4 symmetry. The transformation from a 
stripe SDW state to a state which preserves C 4 symmetry has recently been observed in 
Bai_a;Naa;Fe 2 As 2 near the end of the SDW region^^. 

B. pre-emptive spin-nematic order 

I now analyze a possibility that Z 2 symmetry between X and Y directions gets broken 
before the system develops a stripe SDW order. To analyze this possibility, I include 
fluctuations of the Ax,y fields, introduce the collective Ising-nematic bosonic variable 
cj) cc Ax — Ay together with tjj oc A\ + Ay, integrate over Ax and Ay, and obtain an 
effective action in terms of 0 and 7 /'. I analyze this action and check whether the system 
develops an instability towards (0) 7 ^ 0 before (Ax) or (Ay) becomes non-zero (see Fig. 
21 ). 

That the action (54) can potentially lead to a preemptive Ising-nematic instability is 
evident from the presence of the term g (A^- — Ay)^, which can give rise to an ordered 
state with (Ay) — (Ay) 7 ^ 0 in a way similar to how the term in the Hamiltonian 

(49) gives rise to a state with non-zero {si^) 7 ^ 0. The pre-emptive Ising-nematic insta¬ 
bility, however, does not appear in the mean-held approximation simply because when 
magnetic huctuations are absent, a non-zero (Af) 7 ^ 0 appears simultaneously to (Aj) 7 ^ 0, 
once To changes sign. However, it may well happen once we go beyond mean-held and 
include magnetic huctuations. 

To study a potential preemptive Z 2 transition, I need to introduce collective variables 
of the helds Ax and Ay. Let me introduce auxiliary scalar helds 0 for Ay — Ay and 
0 for Ay -|- Ay- The held ip always has a non-zero expectation value {ip) 0 0, which 
describes Gaussian corrections to the magnetic susceptibility Xi,l Eq. 58. Meanwhile, 
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the field 0 may or may not have a non-zero expectation value. If it does, it generates a 
non-zero value of (A^ — Ay) and the system develops an Ising-nematic order. 

The effective action in terms of 0 and 'ijj is obtained by using again the Hubbard- 
Stratonovich transformation of Eq. (51), but this time the variable x is either ip oc 
A^ + AI- or 0 oc — Ay. Applying this transformation and integrating over fluctuating 
fields Ax and Ay, I obtain the effective action in terms on 0 and ip in the form 

Se» +-'py-<l>P] (57) 

As it is customary for the analysis of fluctuating fields Ax and Ay, we extended the mass 
term tq to include spatial and time variations of Ax,v: 

ro ^ XPi = ro + 7 (58) 


where = 2TiTn is the bosonic Matsubara frequency. 

This action can be straightforwardly analyzed in the saddle-point approximation 
(9S'efr [0,0]/(90 = ^Aefr [0,0]/i90 = 0 (for justification see Ref.®^). Differentiating, I 
obtain two non-linear coupled equations for 0 and ip: 

0 _ f rp + Ip + + j\un\ 

U A (ro + 0 + + 7 |z/„|)2 - 02 


0 ^ f _ 0 _ 

g 4 (ro + 0 + + 7|z/„|)2 _ 02 


(59) 


The full solution of these equations at various temperatures and in different dimensions 
is presented in Ref.®^. The key point is that, for positive g, 0 becomes non-zero at a 
higher temperature (T„) than the one {Tgdw) at which SDW order sets in. In the interval 
Tn > T > Tsdw, (Ax — Ay) becomes non-zero, while (Ax) = (Ay) = 0. Such an order 
breaks C 4 lattice symmetry down to C 2 and is often called Ising-nematic order. 


In Fig. 22 I present the phase diagram for anisotropic 3D system. The transition to an 
Ising-nematic state can be either second-order, or first order. A strong first-order nematic 


transition may instantly induce SDW order. 


C. consequences of the Ising-nematic order 

Because spin-nematic order breaks C 4 lattice rotational symmetry, it couples linearly 
to any other parameter which breaks the same symmetry, such as orbital and structural 
order parameters. Then, once spin-nematic order becomes non-zero, it acts as "external 
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Figure 22: Evolution of the character of the magnetic and nematic transitions in 
the spin-driven nematic theory. The control parameter is the inverse nematic coupling g, 
which changes as function of various control parameters within an itinerant scenario (arrows). 
Second-order (first-order) lines are denotes by solid (dashed) lines. Regions (I)-(III) correspond 
to those of the phase diagram in Fig. 1. The arrows show how the nematic order parameter g 
is expected to evolve with doping, disorder, pressure, and elastic coupling. From [63]. 

field" to the two other parameters and induces non-zero values of both of them. As a 
result, below the fermionic dispersion becomes anisotropic, the occupations of d^z and 
dyz orbitals become non-equal, and also the lattice constants a and b along the x and y 
directions of the Fe-plane, respectively, become non-equal. I refrain to discuss this issue 
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in more detail here and direct a reader to a recent review®^. The development of the 
Ising-nematic order also gives rise to an increase of the magnetic correlation length, what 
in turn gives rise to a pseudo-gap-type behavior of the fermionic spectral function. 

VII. THE STRUCTURE OF THE SUPERCONDUCTING GAP 

I now turn to superconductivity. Like I did for SDW order, I assume that renormal¬ 
izations captured within pRG are already included into consideration and consider an 
effective low-energy model with effective pairing interactions in the band basis. In the 
discussions in this Section I follow Refs. [31,32,84,101,108,112,136,150]. 

A. The structure of s—wave and d—wave gaps in a multi-band SC - general 
reasoning 

In previous sections I assumed that the interactions in the particle-particle channel 
(the dressed G 3 and G 4 terms) are independent on the angles along the hole and electron 
FSs. In this situation, the only option is an s—wave gap, which changes sign between 
the FSs, but is a constant along each FS. Now I consider realistic models in which the 
interactions in the band basis are obtained from the underlying multi-orbital model. These 
interactions generally depend on locations of fermions along the FS. 

I hrst display general arguments on what should be the form of the gap in different 
symmetries and on different FSs. I show that an s—wave gap generally has angle depen¬ 
dence and may even have nodes, while a d-wave gap, which is normally assumed to have 
nodes, may in fact be nodeless on electron FSs. 

A generic low-energy BCS-type model in the band basis is described by 

n = J2ei{k)al^ai]^+ ^ p)ajkaLk«iP«i-P (60) 

i,k 

The quadratic term describes low-energy excitations near hole and electron FSs, labeled 
by i and j, and the four-fermion term describes the scattering of a pair [k —k ),) on 
the FS i to a pair (p —p |) on the FS j. These interactions are either intra-pocket 
interactions (hole-hole Uhihi or electron-electron UeiSi ), or inter-pocket interactions (hole- 
electron Uejhi, hole-hole Uhi^^hj, and electron-electron Ugi^e^)- 

Assume for simplicity that the frequency dependence of F can be neglected and low- 
energy fermions are Fermi-liquid quasiparticles with Fermi velocity v^p- In this situation. 


56 



the gap A(/c) also doesn’t depend on frequency, and to obtain one has to solve the 
eigenfunction/eigenvalue problem; 

\iAi{k) = - J r(kj., Pf)A^(p) ( 61 ) 

where Aj are eigenfunctions and Aj are eigenvalues. The system is unstable towards pairing 
if one or more A* are positive. The corresponding Tc^i scale as Tc^i = Although 

Aj are generally different for different i, the exponential dependence on l/A* implies that, 
most likely, the solution with the largest positive A* emerges hrst and establish the pairing 
state, at least immediately below T^. 

Like I discussed in the Introduction, the pairing interaction U{k,p) can be decomposed 
into representations of the tetragonal space group (one-dimensional representations are 
Aig, Big, B 2 g, and A 2 g). Basis functions from different representations do not mix, but 
each contains inhnite number of components. For example, s—wave pairing corresponds 
to fully symmetric Aig representation, and the s—wave {Aig) component of U{k,p) can 
be quite generally expressed as 

£/<'»)(*:,?) = U.(k,p) (62) 

m,n 

where T^(/c) are the basis functions of the Ai^ symmetry group: 1, cos kx cosky, 
coskx + cosky, etc, and are coefficients. Suppose that k belongs to a hole FS and is 
close to /c = 0. Expanding any wave function with Ai^ symmetry near k = 0, one obtains 
along |k| = kp, 

+ &mCOs40fc COS -f- ... (63) 

where <pk is the angle along the hole FS (which is not necessary a circle). Similarly, for Big 
representation the wave-functions are cos kx — cos ky, cos2kx — cos2ky, etc, and expanding 
them near k = 0 one obtains 

= a*^cos2(j)k + b*^cos6(j)k + c*^cosl0(j)k + ... (64) 

There are no fundamental reasons to expect that bm, Cm or b*^, c*^ are much smaller 
than Om or a*^, but sub-leading terms are often small numerically. I assume that this 
is the case and neglect subleading terms, i.e., assume that s—wave interaction between 
fermions on the hole FSs can be approximated by an angle-independent {k,p) = Uhihj 
{hi label different hole FSs), while d—wave {Big) interaction can be approximated by 

Ut^hi{k,p) = Uhihi cos20fccos20p. 
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The situation changes, however, when I consider the pairing component involving 
fermions from electron FSs. Suppose that k are still near the center of the BZ, but p 
are near one of the two electron FSs, say the one centered at (0,7r). Consider all possible 
'^n{p) with Aig symmetry A simple experimentation with trigonometry shows that there 
are two different subsets of basis functions: 

subset I : 1, cos Px cos Py, cos2px + cos2py... 

subset II ; cospx + cospy, cosSpx + cosSpy... (65) 

For a circular FS centered at (0,7r), the functions from subset I can be again expanded 
in series of cos4/0p with integer 1. The functions from subset II are different ~ they all 
vanish at (0,7r) and are expanded in series of cos(2(;/)p + 4/0p) (the hrst term is cos20p, 
the second is cosG^p, etc). For elliptic FS cosd/^p and cos(20p + d/^p) terms appear in 
both subsets. In either case, the total 

'J'm(p) = dm + &mCOs40p + CmCOs80p + ... 

+am cos 20p + bm cos 60p + Cm COS lO0fc + ••• (66) 

For the other electron FS, T^(p) is the same, but momentum components Px and p^are 
interchanged, hence the sign of all cos(2(;/) + d/^p) components changes. 

Let’s make the same approximation as before and neglect all components with I > 0. 
Then 

'F^(p) =am±am cos 20P (67) 

where the upper sign is for one electron FS and the lower for the other. It is essential 
that the angle-independent term and the cos 2(j)p term have to be treated on equal footing 
because each is the leading term in the corresponding series. Combing (67) with the fact 
that T^(A;) can be approximated by a constant, we obtain a generic form of the s—wave 
component of the interaction between fermions near hole and electron FSs 

(1 + 2ae,h COS 20p^, + ...) 

Ue2,hS^:P) = f/eA (1 - 2ae,feCOs20p^2 + ....) (68) 

where dots stand for cosd^A,, cosd^p, cos60p, etc terms. 

By the same reasoning, s—wave components of inter-pocket and intra-pocket interac- 
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tions between fermions from electron FSs are 


= Ue,e (1 + 2aee (cOS + COS J 
+4/3ee COS 20fc^j COS 2(f)p^^ + ... 

^ 62,62 = Ue,e (1 - 2Q;ee (cOS 20fc^2 + COs20p^2) 

+4/3ee COS 20fc^2 COS 20p^2 + ... 

^61,62 (1 + 2aee (cOS 20^^^ - COs20p^J 

-4/3ee COs20fc^j COs20p^2 + ... 


(69) 


Once the pairing interaction has the form of Eqs. (68) and (69), the gaps along the hole 
FSs are angle-independent (modulo cos40 terms), but the gaps along the two electron 
FSs are of the form 

{k) = Ae ± Ae cos 20fc. (70) 

When Ae is small compared to Ae, the angle dependence is weak, but when |Ae| > |Ae|, 
s—wave gaps have nodes at “accidental” values of 0, which differ between the two electron 
FSs. 

A similar consideration holds for dx‘ 2 _y 2 gap. Within the same approximation of leading 
angular momentum harmonics, we have 


and 


Uer,hS^,P) = Ue,hiCOs2(l)hi{l +ae,hCOs2(l)p^^) + ... 

= ^e,fciCOs 20 /,. (-l + de,fcCOs 20 p^J + ... ( 71 ) 


K 

,eA^:P) = 

= Ue,e (1 

+ 2 Q;ee ( 

cos 20 fc^j 

-f cos 20 p^J 



+40ee 

COS 20fc^j 

cos 20 p^j 

+... 

K 

,e2ik,p) = 

= Ue,e (1 

2 o^ee ( 

cos 20 fc ^2 

+ cos 20 p^J 



+40ee 

COS 20fc^2 

cos 20 p ^2 

+... 

K 

,62 = 

= Ue,e (- 

-1 - 2 aee (cos 20 fc, 

- cos 20 p ^2 



+40ee 

cos 20 fc^j 

cos 20 p ^2 

+... 


The solution of the gap equation then yields the gap in the form 

Aj(^^(/c) = A/iCos20fc 

A^'^) {k) = ± Ae + Ae COS 20 fc. ( 73 ) 
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Along the hole FS, the gap behaves as a conventional d—wave gap with 4 nodes along the 
diagonals. Along electron FSs, the two gaps differ in the sign of the angle-independent 
terms, and have in-phase cos 20 oscillating components. When Ag << Ae the two electron 
gaps are simply “plus” and “minus” gaps, but when Ag > Ae, each has accidental nodes, 
again along different directions on the two electron FSs. 

We see therefore that the geometry of the FSs in FeSCs affects the gap structure in 
quite fundamental way: because electron FSs are centered at the k points which are not 
along BZ diagonals, s—wave gaps on these FSs have cos 20 oscillations which one normally 
would associate with a d—wave symmetry, and d—wave gaps have constant (plus-minus) 
components which one would normally associate with an s—wave symmetry. When these 
“wrong” components are large, the gaps have accidental nodes. These nodes may be 
present or absent for both s—wave and d—wave gaps. 

An s—wave gap with nodes in one of the “exotic” options offered by the electronic 
structure of FeSCs. Another “exotic” option is a d—wave state without nodes. In heavily 
electron-doped FeSCs, hole states are gapped, and only electron FSs remain. The d—wave 
gaps on these two FSs have no nodes if cos 20 oscillation component is smaller than 
a constant term, hence the system will display a behavior typical for a fully gapped 
SC despite that the gap actually has a d-wave symmetry. There are even more exotic 
options offered by the actual three-dimensionality of the electronic structure and/or the 
hybridization of the electron FSs due to interaction via a pnictide/chalcogen. Refs. [93, 
94,133,207]. 

A generic analysis of the eigenvalue/eigenfunction problem, Eq. (61), reduces to the 
set of either four (or Eve) coupled equations in either s-wave or d-wave channels: two (or 
three) A’s are the gaps on the hole FSs, and two other A’s are angle-independent and 
cos 20 components of the gaps on the electron FSs. Accordingly, there are either four or 
five different and A^. 

1. Generic condition for a non-zero Tc- 

Before I analyze specific cases of 4 x 4 and 5x5 gap equations, I consider the issue 
whether in the presence of angular dependence of the interactions its is still required for 
superconductivity that the inter-pocket interaction Uhe must exceed the threshold set by 
intra-pocket hole-hole and electron-electron interactions. Interestingly enough, this may 
no longer be necessary. To illustrate this, consider the case of an s—wave pairing in a 
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four-pocket model and assume for simplification that only one hole pocket is relevant to 
the pairing. Then the eigenvalue problem reduces to the set of three equations for A^, A^, 
and Ae {Ae{k) = Ae-|-AeCos2(;/)fc). Solving the set, we hnd three solutions A| {i = 1,2,3). 
In the absence of cos 20 terms in Tij{k,p), Ag = 0, and Af 2 given by 

~{uhh + 2nee) i \l{Ujih ~ 2Mee)^ + ^^/le 

Ai,2 = 2 

I remind that Uij = UijNQ, where Nq is the density of states. Obviously, Uhe has to exceed 
a threshold, otherwise A^ 2 < 0- Once the angle dependent terms in (68-69) become non¬ 
zero, Ag also becomes non-zero, and its sign depends on the interplay between ate, ctee, 
and /dee- In particular, when < UeeUhh (and, hence, Ai ,2 < 0), Ag is positive or negative 
depending on whether or not A > 0, where 

A AUggUfih ^Clgg /dee^ T ^he d~ 2/dee ScT/jeOiee^ 


When the angle-dependence of the electron-electron interaction can be neglected, i.e., 
ttee = /dee = 0, Ag > 0 uo matter what is the ratio of n^g and UeeUhh- In particular, for 

UhhUee > mL and a/ie << 1, 


> 0 (76) 

^^UhhUee-ul^ 

In other words, for one of s—wave solutions, A^ > 0 even if intra-pocket repulsions are 
the largest. The full solution of the 3x3 set with o/gg = /dee = 0 shows that two A’s 
are repulsive and one is attractive for arbitrary u\^/UeeUhh- When the ratio is small, the 
attractive solution is close to (76), when the ratio is large, the attractive solution is close 
to Ag in (74). I illustrate this in Fig. 23 

There is, however, one essential difference between the cases u\^/UeeUhh > 1 and 
u\e/UeeUhh < 1- In the first case, momentum-dependence of the interaction just modihes 
the “plus-minus” solution which already existed for momentum-independent interaction. 
In this situation, the gap along electron FS gradually acquires some cos 20 variation and 
remains nodeless for small ahe- In the second case, the solution with A > 0 is induced 
by the momentum dependence of the interaction, and the eigenvalue corresponding to Ag 
necessary has Ag > Ag, i.e., s—wave gap has nodes along the electron FS [150]. In other 
words, the pairing occurs for all parameters but whether the gap is nodal or not at small 
ahe depends on the relative strength of intra-pocket and inter-pocket interactions. When 
intra-pocket interaction dominates, the gap “adjusts” and develops strong cos 20 compo¬ 
nent which does not couple to a momentum-independent Uee term and by this effectively 
reduces the strength of electron-electron repulsion. 
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Figure 23: (a) The three eigenvalues in the s— wave channel Af as functions of u\^/{ueeUhh) for 
Oee = /3ee = 0 and Uhe = 0.4. For any u\^/{ueeUhh) , one Af is positive (attractive), other two are 
negative. Positive Af corresponds to s± pairing. At small u\^/{ueeUhh) pairing is induced by 
Uhe and the gap has nodes on electron FSs. At large u\^/{ueeUhh) positive Af exists already at 
oihe = 0, and the gap along electron FS has nodes only if ahe is above the threshold. The circle 
marks the area where positive and negative solutions come close to each other. The splitting 
between the two increases with a^e- (b) The regions of nodeless and nodal gap, depending 
on ahe and u\^/UeeUhh- From Ref. [108]. 


The same reasoning holds for the case of two non-equivalent hole FSs, and for 5-pocket 
models, and also for the d—wave channel. For all cases, the solution with Aj > 0 may 
exist even when intra-pocket interactions are the largest, but in this situation the gaps 
must have accidental nodes. The existence or non-existence of the solution at strong 
intra-pocket repulsion then depends on the complex interplay between the prefactors of 
cos 2d terms in electron-hole and electron-electron pairing vertices, see Eq. (75). 


B. How to extract f/jj(k,p) from the orbital model? 

So far, in our discussion etc, are treated as some phenomenological inputs. To 

obtain the actual values of these parameters, one needs a microscopic model. The most 
commonly considered model for FeSCs is an effective 5-orbital model for Fe atoms with 
local intra-orbital and inter-orbital hopping integrals and intra-orbital and inter-orbital 
density-density (Coulomb) repulsions, Hund-rule exchange, and the pair hopping term. 
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Figure 24: Representative fits of the interactions rjj(ky7’, py?) by LAHA for the 4-pocket model. 
Tij are obtained by converting the Hamiltonian, Eqs. (77), (78) from the orbital to the band 
basis. The symbols represent interactions computed numerically for the 5-band orbital model 
using LDA band structure, the black lines are the fits using Eqs. (68)-(72). The fit is for the 
set U = 1.67, J = J' = 0.21, V = 1.46, and /r = 0.08 (all in eV). A positive /r corresponds to 
electron doping, kyr’ in rjj(ky7’, py?) is selected along y direction on either an electron or a hole 
ES (its location is specified on top of each figure), and pyr is varied along each of ESs. The angle 
(f) is measured relative to kx- 


Hint ^ ( U ^ ^ ^ ( JstSis'Sii-\- ^ ^ ( J ^ ^ 

is i,Syt^s ^ 

where riis = Ui^st + ^i 4 - 

The Hamiltonian Hint can be equivalently re-expressed via spin-independent interac¬ 
tions, as 


H 


int 


= YlUrii^s^nisi + ^Tii 

is i^s^t^s 


J 


■^it 


i,s,tj^s 


t 


^isa ^ita 


h E ^'E 


i,s,t^s 


^isa^isa ^ita ^ita 

(78) 


where U = Vst + Jstl‘^- 

The hopping integrals (36 total) are obtained from the £t to DFT band structure. [131] 
For the interaction parameters, the most common approximation is to assume that H, J 
and J are independent of the orbital indices s and t, as long as s 7 ^ f. The model can be 
also extended to include non-local Fe-Fe interactions via a pnictide [151]. 

The bare parameters in (77) and (78) are inter-related due to local spin-rotation in¬ 
variance [24,60], but that invariance is broken if we view (77) and (78) as an effective 
low-energy model in which the interactions are dressed by the renormalizations coming 
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from fermions with energies of order bandwidth. By this reason, in most stndies U, U, J, 
and J' are treated as independent parameters. 



^hihi 

'^h2h2 

'^h\h2 

rt/iie 

^h\e 

^/i2e 

^/i2e 

^ee 


Pee 

NSF 

0.8 

0.76 

0.78 

0.46 

-0.24 

0.4 

-0.30 

0.77 

0.14 

0.09 

SF 

2.27 

2.13 

2.22 

4.65 

-0.34 

2.29 

-0.22 

3.67 

0.15 

0.04 


^hihi 

^h2h2 

^h\h2 

^h\e 

^h\e 


^h2e 

^ee 

CKee 

he 

NSF 

0.7 

0.66 

-0.68 

-0.25 

-0.58 

0.24 

-0.42 

0.11 

-0.5 

0.25 

SF 

1.50 

1.40 

-1.50 

-3.73 

-0.44 

1.44 

-0.32 

1.03 

-0.49 

-0.02 


Table I: Table for s-wave and d—wave parameters for the same set as in Fig. 24. NSF and 
SF mean the bare interaction without the spin-fluctuation component and the full interaction, 
respectively. 

We now need to convert (77), (78) into the band basis and re-express it in the form 
of Eq. (60). This is done by transforming into the momentnm space, introdncing new, 
hybridized operators, which diagonalize the hopping Hamiltonian, and re-expressing the 
interaction terms in (77) or (78) in terms of these new operators. The end result of this 
procedure is the effective Hamiltonian in the band basis which has the form of Eq. (60) 
with 7/ij(k, p) given by 

U„(kX) = i:“.‘'*(-k)Qr(k)Re|rS(k,k')] 

stpq 

xaP(k')a](-k'), (79) 

where [f/f/(k, k') are linear combinations of U, U, J and J, and af is the matrix element 
connecting the original fermionic operator Cp in the orbital basis with the new fermionic 
operator a* on FS i in the band basis. The matrix elements af contain information which 
orbitals mostly contribute to a particular segment of a particular FS [24,25]. Because 
of this, the interaction 7/ij(k, p) in the band basis generally depends on the angles along 
different FSs and contains components in all representations of the tetragonal D^h group. 

The angle dependence of s—wave and d^ 2 _y 2 vertices agrees by symmetry with Eqs 
(68)-(72). What s a’priori unknown is how well the interactions can be approximated 
by the leading angle harmonics, i.e., whether the terms labeled as dots in (68)-(72) can 
actually be neglected. This issue was analyzed in detail in Ref. [84], and the answer is 
affirmative - the leading anhular harmonic approximation (LAHA) works rather well. In 
Fig.24 I show representative hts for a particular set of parameters and in Table 1, in the 
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lines marked NSF, I show Ueh and other parameters, extracted from the £t (NSF stangs 
for "no spin fluctuations", meaning that this is for the bare interaction, without extra spin- 
fluctuation component (see below)). The results somewhat vary depending on the values 
of U, V, J, J', but in general intra-band interactions in the s—wave channel, Uee and Uhh, 
exceed interband Uhe- This is not surprising because Uee and Uhh are essentially Coulomb 
interactions at small momentum transfers, while Ueh is the interaction at large momentum 
transfer, and it should be smaller on general grounds. Only when V = J = J' = 
the interaction in the band basis becomes independent on the momentum [112], i.e., 
Uee = Uhh = Uhe (tliis was termed “Coulomb avoidance” in Ref. [28]). According to Table 
I, intra-band interactions are also larger in the d-wave channel: UhihiUee > u\.^, although 
the reasons why this is the case are not transparent. 



s-wave solution d-wave solution 



Figure 25: Representative case of small/moderate electron doping, when both hole and electron 
pockets are present. Panel a - the FS, panel b - representative fits of the interactions by LAHA 
(the dots are RPA results, the lines are LAHA expressions, Eqs (68)-(72)). Panels c and d - the 
eigenfunctions in s—wave and d—wave channels for the largest A® and X^. From Ref. [84]. 
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Figure 26: The fits of the RPA interactions by LAHA and the structure of s—wave and d—wave 
gaps for the case of heavy electron doping, when only electron FSs are present. From Ref. [84]. 

C. Doping dependence of the couplings, examples 

I now present the results for the gap structure at various doping, obtained within 
LAHA, but including dressing of interactions by spin-fluctuations in RPA (lines marked 
"SF” in the Tables). The results for hole and electron doping differ, and I present them 
separately. I will follow Refs. [31,84]. 


1. Electron doping 

For small and moderate electron dopings, the FS consists of 4 pockets - two hole FS 
at (0, 0) and two electron FSs at (0, vr) and (vr, 0). Typical fits by LAHA, the parameters 
extracted from the hts, and the solutions in s-wave and d-wave channels are shown in 
Fig. 25 and in Table H. It turns out [84] that some system properties are sensitive 
to the choice of the parameters, but some are quite universal. The parameter-sensitive 
properties are the presence or absence of accidental nodes in the s-wave gap (although 
for most of parameters the gap does have nodes, as in Fig. 25) and the gap symmetry 
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itself, because for most of input parameters and dopings and A'^ remain comparable 
as long as both hole and electron FSs are present (see Table II). That d—wave state is 
a strong competitor in 4-pocket systems has been hrst emphasized in Refs. [24,60]. The 
authors of [24] hinted that different FeSCs may have different symmetry even for the 
same topology of the FS. 

The universal observation is that the driving force for attraction in both s-wave and 
d-wave channels is strong inter-pocket electron-hole interaction ^'^d Uhie terms) no 
matter how small the hole or electron pockets are. The gap structure actually changes 
only little with doping as long as both hole and electron pockets are present. 


Table II: Some of the LAHA parameters extracted from the LAHA fit in Figs. (25) and (26) 
for electron doping. Blocks (i) corresponds to Fig. (25), block (ii) corresponds to Fig. (26) (no 
hole pockets). From Ref [84]. 


(i) (h) 


s-wave 


^hie 

^h\e 

'^ee 

^ee 

A. 

^ee 

Q^ee 

A. 


0.8 

0.79 

-0.19 

0.91 

0.05 

0.25 

3.65 

0.20 

0.1 

d-wave 

hh\hi 

^hie 

/lie 

^ee 

^ee 

Ad 

dee 

dee 

Ad 


0.50 

-0.39 

-0.46 

-0.04 

1.5 

0.37 

-2.57 

0.29 

5.9 


Extreme electron doping 

The situation changes qualitatively once the hole pockets disappear (Fig. 26). It is 
clear from Table II that now the d-wave channel becomes the dominant one. Comparing 
the LAHA parameters for the two dopings, we see the reason; once the hole pockets 
disappear, a direct d-wave electron-electron interaction Uee becomes strong and attractive. 
The argument why this happens is as follows: [84] Uge is an antisymmetric combination of 
intra-pocket and inter-pocket electron-electron interactions Fee = '*^fi)rra~'*^fnter- Both nf^ter 
and positive (repulsive), but the sign of Fee depends on the interplay between 

and nf()ra- BSs are present, SF are peaked near q = (0, tt) and (tt, 0), 

which are an equal distance from the relevant momenta q = 0 for nf^tra 
for In this situation, nf,()ra wfnter remain close in magnitude, and Fee is small. 

Once the hole pockets disappear, the peak in the RPA spin susceptibility shifts towards 
(tt, tt) [188] and increases more due to the SF component than ^ negative 
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■^^fntra ~ '^tnter then gives rise to a “plus-minus” gap on the two electron FSs. The gap 
changes sign under —)■ ky and therefore has dx^_y 2 symmetry. This pairing mechanism 
is essentially identical to spin-fluctuation scenario for d-wave pairing in the cuprates [100]. 

There are other proposals for the gap structure at extreme electron doping. The authors 
of Refs. [132,208] argued that the gap symmetry may be nodeless s—wave (equal sign of 
the gap on the pockets at (0, n) and n, 0), if one uses for electron-electron interaction the 
weak coupling version of the Ji — J 2 model. Another proposal for strongly electron-doped 
FeSCs is s’*"’*' pairing driven by orbital fluctuations [149]. And yet another proposaT^^^^^^ 
is that the pairing state in FeSCs with only electron pocket present is s’*"", with the sign 
change between the hybridized electron pockets. Such a state emerges if one includes 
into consideration the hybridization of the two electron pockets. In this novel s’*"” state, 
all electron states are gapped, yet because of sign change, there is a spin resonance at 
momenta which is roughly a distance between the electron pockets in the unfolded BZ®^. 


Table III: Some of LAHA parameters extracted from the fits in Figs. 27 and 28 for hole doping. 
Block (i) corresponds to Fig. 27 (hole and electron pockets are present), block (ii) corresponds 
to Fig. 28 ( no electron pockets). 


(i) (hi) 


s - wave Uhihi 

^hie 

^hie 

^ee 

A. 


'^h\h2 


Ufishs, 

As 

0.0.86 

0.92 

-0.18 

1.00 

0.58 

0.67 

0.8 

0.29 

1.37 

0.13 

d - wave Uhihi 



'Uee 

Ad 


^h\h2 


Uh^hz 

Ad 

0.51 

-0.45 

-0.48 

0.07 

0.31 

0.36 

-0.5 

-0.02 

-0.17 

0.11 


2. Hole doping 

For small and moderate hole doping, the FS contains 5 pockets -two hole pockets at 
(0,0), two electron pockets at (0,7r) and (7r,0), and one more hole pocket at {71,71). Rep¬ 
resentative FSs for hole doping, typical hts by LAHA, the parameters extracted from the 
£t, and the solutions in s-wave and d-wave channels are shown in Fig. 27 and in Table III. 
Just like for electron doping, there are universal and parameter-sensitive features. The 
parameter-sensitive property is again the presence or absence of accidental nodes in the 
s-wave gap along the electron FSs, although for most of the parameters, the gap does 
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Figure 27: Representative case of small/moderate hole doping, when both hole and electron 
pockets are present. Panel a - the FS, panel b - representative fits of the interactions by LAHA 
(the dots are RPA results, the lines are LAHA expressions, Eqs (68)-(72)). Panels c and d - the 
eigenfunctions in s—wave and d—wave channels for the largest A® and A'^. From Ref. [84]. 

not have nodes (see Fig. 27) because the total Uhe increases once it acquires an additional 
contribution Uh^e- 

There are two universal features. First, the s-wave eigenvalue is enhanced relative to a 
d—wave one and becomes the leading instability as long as both hole and electron pockets 
are present. Second, the driving force for the attraction in both s- and d- channels is 
again strong inter-pocket electron-hole interaction {uhe and Uhe terms), no matter how 
small electron pockets are. 


Extreme hole doping 

The situation again changes rapidly once electron pockets disappear, see Fig. 28. Now 
electron-hole interaction becomes irrelevant, and the attractive pairing interaction may 
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Figure 28: The fits of the RPA interactions by LAHA and the structure of s—wave and d—wave 
gaps in for strong hole doping = —O.SOeR), when only hole FSs are present. From Ref. [84]. 


only be due to intra and inter-pocket interactions involving hole pockets. LAHA analysis 
shows [84,136] that, at least for in some range of parameters, there is an attraction in both 
s—wave and d—wave channels, and furthermore Xd ~ A^, see Fig. 28 The near-equivalence 
of s—wave and d—wave eigenvalues was also found in recent unrestrictive RPA study [153]. 
Within LAHA, the attractive is due to strong intra-pocket interaction between the two 
hole pockets centered at (0,0). The s—wave gap then changes sign between these two hole 
pockets. The gap along (vr, tt) pocket is induced by a weaker inter-pocket interaction and 
is much smaller. LAHA neglects cos4n0 gap variations along the hole FSs (i.e., s—wave 
gaps are treated as angle-independent), but the theory can indeed be extended to include 
these terms. The attractive Xd emerges by two reasons. First, the d-wave intra-pocket 
interaction Uh 3 h 3 becomes negative, second, the inter-pocket interaction between the 
two pockets at (0, 0) becomes larger in magnitude than repulsive Uh^hi and Uh 2 h 2 (see 
Table HI). The solutions with A^ > 0 then exist separately for FSs hi ^2 and for h^, the 
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residual inter-pocket interaction just sets the relative magnitudes and phases between the 
(larger) gap at and (smaller) gaps at hi 2 . The d-wave gap with the same structure 
has been obtain in the fRG analysis at large hole doping [152]. 

D. LiFeAs 



Figure 29: Fermi surface of LiFeAs as deduced from the ARPES experiments: (a) shows the 
three-dimensional version of the Fermi surface and (b) and (c) refer to the two-dimensional cuts 
at 1 ^ 2 ! = TT (left) and = 0, respectively. Hole pockets are located a (0,0) and ( 7 r, 7 r) and 
electron pockets are at (± 7 r, 0) and ( 0 ,± 7 r). In case of kz = 0, the two tiny hole pockets hi ^2 
vanish just below the FS and only /is and 61^2 remain, (d) shows the zoomed region of the first 
BZ around the P—point of the BZ with tiny a hole pockets. 

There is a possibility to obtain a more complex behavior even in systems which contain 
both hole and electron pockets. One such example is LiFeAs. Its electronic structure 
contains three hole and two electron pockets, however two T—centered hole pockets have 
strong 3D dispersion and exist only near = vr (see Fig.29 As a result the FS in the 
cross-sections at small kz consists of one hole and two electron pockets (hole 7 pocket and 
electron /3 pockets), while in the cross-section at kz near tt the FS consists of three hole 
pockets and two electron pockets. The orbital content of the FSs for the two a FSs is 
very different from that for other three FSs. Namely, the two F—centered hole pockets {a 
pockets) are made chiefly of d^z and dyz orbitals. The other three FSs are made primarily 
of dxy orbital, with rather small admixture of dxz and dyz orbitals (Refs.^^^’^^®. These 
features indicate that the low-energy electronic structure of LiFeAs consists of two very 
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different subsets. One is made out of quasi-2D 7 and /3 pockets with primarily d^y orbital 
content, and the other is made out of a pockets, which are highly anisotropic along kz 
and are made primarily out of dxz and dyz orbitals. 


s- 

-wave 

^hihi 

F/i2/l2 


Uh\h2 

Uhihs 

F/i2/l3 

F/ije 

^hie 

Uh2e 


CO 

a" 

^hse 

Uee 


Idee 

J: 

= O.OF 

0.92 

0.99 

1.21 

0.95 

0.29 

0.23 

0.28 

-0.34 

0.22 

-0.49 

1.20 

-0.12 

1.20 

-0.12 

0.03 

J: 

= 0.H7 

0.99 

1.09 

1.14 

1.03 

0.16 

0.10 

0.15 

-0.85 

0.08 

-1.65 

1.14 

-0.13 

1.14 

-0.13 

0.04 

J: 

= 0.31/ 

1.14 

1.28 

1.02 

1.20 

-0.09 

-0.15 

-0.12 

1.58 

-0.18 

1.10 

1.02 

-0.16 

1.03 

-0.17 

0.05 


Table IV: LAHA projected interactions in the s—wave channel for = vr. The energies are in 
units of U. 

S-Wave O/ig/jg Uh^e Oih^e Uee OLee f^ee 

J = O.OU 1.53 1.36 -0.14 1.22 -0.13 0.04 

J = 0.1U 1.48 1.30 -0.16 1.16 -0.14 0.04 

J = 0.3U 1.40 1.19 -0.20 1.04 -0.18 0.06 

Table V: LAHA projected interactions in the s—wave channel for kz = 0. The energies are in 
units of U. 

The results^^® for the interactions in s—wave channels within LAHA are shown in 
Tables IV and V. For the model described by Eq. (78) with U' = U — 2J and J' = J, we 
clearly see that the two subsets are nearly separated for all J/U. This near-separation 
opens up a novel possibility for the structure of s’*"” gap. Namely, superconducting gaps 
consistent with the structure of interactions in Table IV are 

Afti(0) = Ahi 

^h2{4>) = A/12 

A. 3 ( 0 )=A ,3 (80) 

Aei {0) = Ae + Ae COS 26 
Aea (^) = Ae — Ae COS 26 

In a "conventional" s’*"” gap structure, the gaps on the three hole pockets are of the same 
sign. Here, superconductivity within the subset of the two a pockets is primary due to 
inter-pocket repulsion between fermions near these pockets. When this repulsion exceeds 
inter-pocket repulsion, it gives rise to sign-changing s’*"" superconductivity between these 
two pockets. In Fig. 30 I show the gap structure obtained for the parameters from Table 
IV and how it evolves when I artihcially increase the interaction between a and /d pockets. 
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When inter-subset interaction is strong, I obtain a conventional s’*"” superconductivity, 
with the same sign of the gap on all three hole pockets^^^. However, for small/moderate 
coupling between the a and /3 — 7 subsets, we see from Fig. 30 that the gaps on the two 
a pockets are of opposite sign. 

Another novel structure of an s—wave gap (termed as "orbital antiphase state") has 
been suggested in Ref.^^®. In this state, superconducting gap has the same sign on the two 
a pockets but changes sign between a pockets and 7 pocket. Such a state occurs if the 
coupling between the two subsets is strong and predominantly involves a — 7 interaction. 

E. Superconductivity which breaks time-reversal symmetry 

Several groups argued recently that multi-orbital character of FeSCs is an ideal play¬ 
ground to search for a truly novel spin-singlet superconductivity which breaks time- 
reversal symmetry (TRS). Spin-triplet superconductivity with broken time-reversal sym¬ 
metry (pa, state) has likely been found in Sr 2 RuO^^^, which represents a solid-state 
analog of superfluid 3fjg54,i96^ spin-singlet d + id state has not yet been observed 

experimentally. Such a state was once proposed as a candidate state for high Tc cuprate 
superconductors^®^, but later gave way to a more-conventional TRS-preserving d—wave 
state. A TRS breaking d-\-id superconductivity has been recently predicted for fermions 
on a hexagonal lattice (e.g., graphene) near van-Hove doping^^^’^^®. 

For FeSCs, one proposal is to explore the region where s—wave and d—wave pairing 
channels are competitive in strength, and there is a transition from one pairing symmetry 
to the other, as one varies the parameters. In the intermediate regime, the system very 
likely falls into an intermediate s + id state, with a broken 

There is an even more exotic possibility to get a broken TRS state in an s-wave 
superconductor^^"^^^"^^. Consider as an example a system at extreme hole doping, like 
Bai_a;Ka,Fe 2 As 2 at a; ~ 1 and assume that the superconducting order is s—wave, with the 
sign change of the gap between the two F—centered hole pockets. Once hole doping gets 
smaller and electron pockets appear, the system eventually develops a "conventional" s’*"” 
superconductivity in which the gaps on the two F—centered hole pockets have the same 
sign. According to theory^^^, the system evolution with decreasing x from a novel to a 
conventional s^~ order may go through an intermediate state in which the relative phase 
(j) between the gaps on the F-centered hole pockets gradually evolves from 0 = ± 7 r in the 
novel s-wave state to 0 = 0 in a conventional s^~ state (Fig. 31). In between, the system 
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Figure 30: Evolution of the gaps (a) and the largest eigenvalue (b) for the s—wave solution, 
with the coupling between the subset A (the two hole a-pockets) and the subset B (the outer 
hole 7 -pocket and the two electron /3-pockets). From Ref. [138] 

selects either 0 or —0, which are related by time-reversal transformation, i.e. it breaks 
time-reversal symmetry (an s ± /s state). 

An intermediate state with broken time-reversal symmetry is also expected in 
Bai_a;Ka,Fe 2 As 2 at x < 1 if the snpercondncting state in KFe 2 As 2 is d—wave. In this 
sitnation, the system must transform from a d—wave at x = 1 to an s-wave at a smaller 
X, and, like I just said, this normally involves an intermediate s±/d phase. Another recent 
proposal for TRS broken superconducting state in FeSCs is s + it superconductivity in 
the co-existence phase with SDW (see Ref. [218]). 


VIII. EXPERIMENTAL SITUATION ON SUPERCONDUCTIVITY IN FESCS 

As of today, there is no “smoking gun” experiment which would carry the same weight 
as phase-sensitive measurements of da, 2 _y 2 gap symmetry in the cuprates [58]. Still, there 
is enough experimental data to minimize the number of possible gap structures. 

As we discussed in the previous section, theoretically proposed gap symmetry and 
structure can be different for weakly/moderately doped systems with hole and electron 
FSs and for strongly doped systems where FSs of only one type are present. It is then 
instructive to consider weak/moderate and strong doping separately. 
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A. Moderate doping, gap symmetry 


The candidates are s-wave (either or or gap. The two behave very 

differently along the hole FSs centered at (0,0) - s-wave gap is nodeless with cos40 
variations, while d-wave gap has nodes along kx = i/cy. ARPES measurements, both from 
synchrotron [88,154-156] and using laser light [66], show quite convincingly that the gap 
along hole FSs is nodeless in both hole and electron-doped FeSCs. This unambiguously 
selects an s-wave. Additional evidence in support of s-wave pairing comes from very flat 
low-T behavior of the penetration depth in the highest Tc 1111 FeSCs systems [157]. 

B. Moderate doping, vs 

The distinction between and gaps is a more subtle issue, particularly given that 
both belong to the same Aig representation and also because in general Aig gap on electron 
pockets may have strong oscillating component. In general, the gaps on electron and hole 
FSs have non-equal magnitudes, and the issue whether the gap is or s’*""'' reduces to 
whether the gap averaged over an electron FS has the same sign or opposite sign than the 
gap averaged over a hole FS. This is not a fundamental symmetry issue and, moreover, 
when cos 20 oscillations are strong, one may switch from equal to opposite signs of the 
averaged gaps by a small change of parameters [108] or by adding impurities. [158] Still, 
when oscillations are not very strong, whether the eigenfunction has or s’*"’*' character 
is essential because it determines, to a large extent, whether the pairing is driven by spin 
or by orbital fluctuations (see Sec.VI). 

The experimental data most frequently cited in support of gap is the observation 
of a magnetic resonance in neutron scattering [15,159]. If, as many researchers believe, 
the resonance is a spin exciton, it exists at a momentum Q if the gaps at FS momenta 
kp and kp + Q are of opposite sign. Experimentally, in most FeSCs the resonance is 
observed [15,159] near Q = (7r,7r) in the folded BZ, which in this zone is precisely the 
distance between electron and hole FSs. The excitonic resonance then exists if the gap 
changes sign between hole and electron pockets and does not exist if the gap doesn’t 
change sign. A similar reasoning has been used in identifying the the resonance seen in 
the cuprates with a fingerprint of dx^-y^ gap symmetry [160] 

The neutron peak is the resonance if it is narrow and is located below twice the gap 
value. The argument made by the supporters of s’*"’*' scenario [61] is that the observed 
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neutron peak is more broad than the resonance seen in the cuprates, and that there is 
no hrm evidence that the peak energy is below 2A for the minimum gap. For gap 
structure, there is no resonance, but there is a redistribution of the neutron spectral weight 
immediately above 2A what gives rise to a local maximum in the magnetic structure 
factor [61,161,162], Still, the majority of researchers do believe that the observed neutron 
peak is a resonance, and the fact that it is quite broad is at least partly due to cos 2(f) gap 
variations along the electron FSs [162], 

Another rather strong evidence in support of gap is the observed variation of the 
quasiparticle interference pattern in a magnetic held [67] although the interpretation of 
the data has been subject of debates [163]. It was also argued [164] that the very presence 
of the co-existence region between SC and stripe magnetism in FeSCs is a hngerprint of 
an gap, because for s’*"’*' gap a hrst order transition between a pure magnetic and a 
pure SC state is a much more likely scenario. 

C. Moderate doping, nodal vs no-nodal gap 

Let’s assume that the pairing is driven by spin huctuations and the gap has struc¬ 
ture. In 2D scenario, such gap has cos 20 variations along electron FSs, which, according 
to theory, can be rather strong, particularly in electron-doped FeSCs. Experimental data 
show that, whether or not the gap is nodeless or has nodes, depends on the material, on 
the doping, and on whether SC co-exists with SDW order. 

1. Hole doping 

For hole-doped FeSCs (e.g. for Bai_a;K 2 ;Fe 2 As 2 ) the data indicate that the gap is 
nodeless, away from the co-existence region. This is consistent with the theory (see Sec. 
VII). ARPES experiments do not show any angular variation of the gap along both hole 
and electron FSs [88,154], but it is not entirely clear whether ARPES can at present 
distinguish between the gaps on the two electron FSs which in folded zone are both 
centered at (tt, tt). Thermal conductivity data show that k/T tends to zero in the limit of 
T = 0, in line with what one should expect for a nodeless SC [165]. Specihc heat data also 
show non-nodal behavior [166]. The interpretation of the penetration depth data requires 
more care as the data do show a power-law behavior A(T) — A(0) oc T“ with a ~ 2 (Refs. 
[167]). Such a behavior is expected for a SC with point nodes, but it is also expected in a 
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wide range of T for a nodeless SC in the presence of modest inter-band scattering by 
non-magnetic impurities [168]. Penetration depth measurements on artihcially irradiated 
samples [169] support the idea that the gap is nodeless and power-law T“ behavior of 
A(T) — A(0) is due to impurities. 

2. Electron doping 

For electron-doped FeSCs, e.g., 122 materials like Ba(Fei_a:Coa;) 2 As 2 or 1111 materials 
like NdFeAsOi-xFj;, ARPES shows no-nodal gap along hole FS [155,156], but there are 
no data on the gap along each of the two electron FSs. At optimal doping, the data on 
both thermal conductivity [170,173] and penetration depth [173,174] are consistent with 
no-nodal gap However, the data for overdoped Ba(Fei_j,Coa;) 2 As 2 indicate that gap nodes 
may develop; the behavior of A(T) becomes more steep, and k/T now tends to a hnite 
value [170], expected for a SC with line nodes. The data also show \/H behavior of n in 
a magnetic held [170] expected for a SC with line nodes [171], but it was argued that the 
behavior resembling can be obtained even if gap has no nodes [172]. There is also 
clear anisotropy between in-plane conductivity and conductivity along direction, what 
was interpreted [170] as an indication that the nodes may be located near particular k^- 
Specihc heat data in overdoped Ba(Feo. 9 Coo.i) 2 As 2 were also interpreted as evidence for 
the nodes. [175] 

The development of the nodes in gap upon electron doping is in line with the theory. 
The farther the system moves away from the SDW phase, the weaker is the increase of 
intra-band electron-hole interaction and hence the stronger is the competition from intra¬ 
band repulsion. As I discussed in Sec.Ill), the gap adjusts to this change by increasing 
its cos 20 component in order to effectively reduce the effect of the intra-band repulsion 
in the gap equation. 

There is also experimental evidence for cos 20 gap oscillations from the observed oscilla¬ 
tions [176] of the held-induced component of the specific heat C{H, T) in superconducting 
FeTei_ 2 ;Sea; {x ~ 0.5). The measured C{H,T) oscillates with the direction of the applied 
held as cos 40. In theory, such an oscillation is related to the behavior of A^(0) (Ref. 
[193]), hence cos20 gap oscillations in A lead to cos40 oscillations in C{H,T). The ob¬ 
served held and temperature dependence of the prefactor for cos 40 term are consistent 
with the idea that the oscillations are caused by cos 20 term in A. These data were also 
interpreted as evidence for no-nodal gap because if cos 20 gap oscillations were strong and 


77 



the gap had nodes at accidental points, the behavior of would be more complex than 
the observed a + b cos 40. 

For LiFeAs, which is undoped but has FS structure similar to electron-doped FeSCs, 
no-nodal behavior has been observed in ARPES [177], specihc heat [178], penetration 
depth [179] and NMR [90,180] measurements. An a-|-6 cos 40 variation of the gap on the 
hole 7 pocket, consistent with s—wave superconductivity, has been observed^^^, together 
with c + d\ cos 201 variations of the gap on electron pockets^^®. The latter is precisely what 
is expected theoretically for an s-wave superconductor when the hybridization between 
the two 0 pockets is weak. The gap on the a pocket probably also has angle dependence, 
but the pocket is too small to detect it in ARPES measurements. 

3. Co-existence region with SDW 

Taken at a face value, thermal conductivity and penetration depth data indicate that 
the gap becomes nodal deep in the co-existence regime in both hole-doped and in electron- 
doped FeSCs. The most striking evidence comes from thermal conductivity [165,170] - 
in the co-existence regime k/T tends to a hnite value at T —)• 0 and shows y/H behavior, 
both typical for a SC with line nodes. From theory perspective, the gap remains nodeless 
near the onset of co-existence with SDW if it was nodeless outside f the co-existence 
phase^®^, however, deep in the co-existence phase angular variation of the gap increase 
due to FS reconstruction (Ref. [219]). 

4- Isovalent doping 

Electron or hole doping is not the only way to change the properties of FeSCs. Another 
route is to replace one pnictide with the other. The most common replacement is As —)• P. 
P-containing materials include the very Erst FeSC - LaFeOP, with < 5K (Ref. [182]), 
the family BaFe 2 Asi_j;Pa; with the highest Tc around 30K (Ref. [183]), and LiFeP [184]. 
Penetration depth, thermal conductivity, specific heat, and NMR data [185] in these 
materials all show the behavior consistent with line nodes. In particular, k scales linearly 
with T at low T and displays \/H behavior in a magnetic held, and A(T) — A(0) is also 
linear in T down to very low T. Laser ARPES data show [66] that the gap along FS is 
nodeless, so the nodes likely are located on electron FSs. 

On general grounds, the existence of the nodes on electron FSs is in line with theory 
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Figure 31: Qualitative phase diagram for strongly hole-doped FeSC. I model the doping de¬ 
pendence by varying the ratio of inter-pocket electron-hole and hole-hole interactions Uhe/uhh- 
The -|— state has gaps of opposite signs on the two GCP’s and no gap on electron pockets, the 
state is an ordinary s± state, in which the gaps have opposite signs on hole and electron 
pockets. The state with broken TRS is in between the two. The gap structures are pictorially 
presented inside each region by vectors placed inside the circles. The magnitudes of the vectors 
represent |Aj| and the angles represent the phases. Cases (a) and (6) are for equal and non-equal 
intra-pocket interactions and Uh 2 ) for the two hole pockets, respectively . For (a), the state 
with broken TRS starts right at Tc and extends into a finite range at T = 0. For (6), this state 
splits off from the Tc line and is only accessible at lower temperatures, while immediately below 
Tc the -|— state gradually evolves into the state as Uhe/u^h increases. 

predictions particularly as BaFe 2 Asi-^jP^, has the same structure of 4 cylindrical FSs as 
electron-doped FeSCs for which nodes are most likely. It has been argued [60] that a 
replacement of As by P changes the hight of a pnictide with respect to Fe plane, what 
effectively reduces inter-pocket electron-hole interaction, in which case the gap develops 
nodes to reduce the effect of intra-pocket repulsion. However, this argument is only 
suggestive, and it is not entirely clear at the moment why all P-based FeSCs have nodes. 
One way to analyze this semi-quantitatively is to study the correlation between 2A/Tc on 
the hole FS and the presence of the nodes on electron FSs. This study shows [186] that 
from this perspective P-based FeSCs are indeed the “best case” for the gap nodes. 

Another open issue is the location of the nodes along z- direction. Oscillations of ther¬ 
mal conductivity with the direction of a magnetic held have been measured recently [187], 
and cos 40 component of these oscillations has been interpreted using the modihed 2D 
form of the gap on an electron pocket Ae(/c^) = Ao(l + (y{kz) cos20). The best ht to the 
data yields a{kz) > 1 for some kz and a{kz) < 1 for others, in which case the nodes form 
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patches along k^- This gap structure has been reproduced in microscopic calculations^®^, 
but whether this is the only explanation of the data is unclear. 

It is still possible, though, that the nodes are located on a hole FS, near particular kz, 
as some of 3D theories suggest [94], Another possibility, which is also not entirely ruled 
out, is that the system behavior near the surface, probed by ARPES, is not the same 
as in the bulk. The probability that this is the case is not high, though, because some 
ARPES data have been obtained using a laser light which probes states located farther 
from the surface than in conventional synchrotron-based ARPES. 

D. Strongly doped FeSCs 

1. Electron doping 

Strongly electron doped materials are represented by a family of A 2 ;Fe 2 -j;Se 2 (A = 
K,Rb,Cs) [10,11] of which K 0 . 8 Fe 1 . 7 Se 2 is the most studied material. Tc in Aa;Fe 2 -j,Se 2 
is rather high, almost 40K. ARPES shows [11] that only electron FSs are present in 
Aj,Fe 2 -ySe 2 , while hole pockets are at least 60meR from the FS, although hole dispersion 
above 60meK is still clearly visible in ARPES. Two electron FSs are at ( 0 , 7 r) and ( 7 r, 0 ), 
like in other FeSCs, and there is, possibly, another electron FS at (0,0). RPA, LAHA 
and fRG calculations for these systems that the gap should have 

a d-wave symmetry, at least for the case when the FSs are only at ( 0 , 7 r) and ( 7 r, 0 ). A 
d-wave symmetry in this situation means that the gaps on the two electron FSs behave as 
Ao(±l-|-a cos 20), and all calculations yield a < 1, i.e., no nodes (neglecting 3D effects). 
One theoretical alternative is s’*"’*' symmetry by one reason [149] or the other [132,208], 
another is s’*"" state between electron pockets®®’^®®’^®^. At present, both ARPES®® and 
specihc heat data^®® point that the gap is nodeless, at least for most of kz values. Of 
particular relevance here are ARPES data on a small electron pocket centered at = tt 
and kx = ky = 0. These data show^®^ that the gap has no nodes, and, taken at a face 
value, rule out d—wave. At the same time, neutron data clearly show®® spin resonance, 
which, if interpreted as spin exciton®®®, requires a sign change of the gap. Both ARPRS 
and neutron data and recent Raman data^^® are consistent with the novel s®"” gap, but 
more studies are needed to verify whether this state is the right one for Aa;Fe 2 -j,Se 2 
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2. Hole doping 


The case of extreme hole doping is represented by KFe 2 As 2 (T^ = 3K), which is at 
the opposite end from parent BaFe 2 As 2 in the family of Ka;Bai_a;Fe 2 As 2 . According to 
ARPES [86], this system has no electron pockets. It contains hole pockets at (0,0) and 
additional hole pockets around (7r,7r), but whether the latter play any role for supercon¬ 
ductivity is not clear at the moment. 

Both thermal conductivity and penetration depth measurements clearly point to nodal 
behavior [89]. There is, however, no “smoking gun” symmetry-sensitive measurement, so 
whether the gap is a d-wave or an s-wave with nodes due to strong cos 40 gap component 
on one of the FSs remains an open issue. Recent study of Tc under pressure has found^°^ 
that Tc initially decreases with pressure initially, and then suddenly changes trend above a 
critical pressure Pc and start increasing. This is a strong indication of the near-degeneracy 
between different pairing states in KFe 2 As 2 . These states can be s and d, or different 
s—wave states. Like I said before, some theorists suggested mixed states, like s+id or s-fis, 
which break TRS. Zero-held /xSR measurements so far have not detected spontaneous 
internal magnetic helds, expected for s + id state^^^. Whether such helds are generated 
for s + is state is less clear. This is an active held of research and I refrain from discussing 
it in more detail. 

3. FeTei-xScx 

There has been high interest recently in the properties of Fe-chalcogenide FeTei_j,Sea;. 
The parent compound FeTca, is a SDW metal, but with high magnetic moment and 
magnetic order diherent from that in Fe-pnictides (see^^^ for details). Superconductivity 
emerges around x = 0.5, and magnetic huctuations and FS-geometry at these x do not 
diher substantially from Fe-pnictides, and it is reasonable to expect that magnetic huctu¬ 
ations may mediate s’*"” supeconductivity. This, however, has not been studied in detail 
yet. On the other end, at a; = 1, the system behavior is quite unusual - the structural tran¬ 
sition sets in at around 80K (see^^^), well before magnetic huctuations develop. Structural 
order increases as T decreases, but changes trend below superconducting Tc ~ 9, what 
clearly shows that structural order competes with superconductivity. The existence of 
structural transition without strong magnetic huctuations fueled speculations that struc¬ 
tural order in FeSe may rehect spontaneous orbital order, i.e., orbital huctuations are "in 
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the driver’s seat". If this is the case, one should expect superconductivity in this 
material. Experimental studies in the superconducting state of FeTei_a:Sea; are called for 
to resolve this issue. There is also an interesting and highly unusual system behavior 
under pressure - structural transition temperature goes down and magnetic fluctuations 
rapidly develop^^^. 

A highly unusual behavior has been detected in thin Elms of FeSe. ARPES measure¬ 
ments detected only electron pockets^^®, like in Ka,Fe 2 -j;Se 2 . Like in Ka,Fe 2 -j;Se 2 , is 
rather high, ~ 60K (even higher T^. have been reported^^^). Whether FeSe (and, more 
generally, FeTei-^^Sea; for x > 0.5) falls into the same category as Fe-pnictides remains to 
be seen. 

E. Summary 

Overall, the agreement between itinerant theory and experiment with respect to the 
type of SDW order, the interplay between the nematic order and magnetism (e.g., the 
normal state phase diagram as a function of doping) and the symmetry and structure of 
the superconducting gap is reasonably good. With respect to superconductivity, itinerant 
approach predicts that the gap is s’*"” in most systems, with angular variation, chiefly 
on electron pockets, and with accidental nodes in some systems. A more complex gap 
structure emerges in systems with strong hole or strong electron doping. 


IX. CONCLUSION 

he analysis of the normal state behavior and superconductivity in FeSCs is a fascinating 
subject because of multi-orbital/multi-band nature of these materials. This review is 
an attempt to present a coherent picture of itinerant scenario for FeSCs. I discussed 
the SDW magnetism, which in most FeSCs corresponds to stripe order, the pre-emptive 
nematic phase, and the origin of superconductivity and the symmetry and structure of 
the superconducting gap in different classes of FeSCs and at different doping levels. It 
is safe to say that there is no major disagreement between theory predictions within the 
itinerant scenario and the experiments. This by no means implies that FeSCs fall into a 
class of "weakly coupled Fermi liquids". There are numerous indications that the coupling 
is strong and is not that much different from that in the cuprates, where, we know, Mott 
state does develop near half-fflling. Nevertheless, most of parent compounds of FeSCs 


82 



are metals, and, in my view, that low-energy physics of FeSCs is adequately captured 
within a moderate coupling itinerant scenario. Up to what temperature/energy one can 
extent the itinerant approach is another question. This scale varies from material to 
material, but, still, is larger than the scales associated with superconducting Tc and is 
likely larger that SDW transition temperature Tgdw and the temperature Tn associated 
with the development of a nematic order. The physics at higher temperature/energies is 
outside the validity of the itinerant approach. 
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